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PREFACE 


The  problem  of  guiding  a  submerged  submarine  requires  the  use  of  es¬ 
sentially  new  navigation  techniques.  One  such  method  involves  the  applica¬ 
tion  of  completely  self-contained  inertial  systems  such  as  are  being  developed 
for  long-range  all-weather  aircraft  operations.  Accordingly,  Instrumentation 
Laboratory  of  Massachusetts  Institute  of  Technology  has  undertaken  the  study 
of  the  application  of  inertial  techniques  to  the  navigation  of  submarines. 

This  study  is  conducted  under  contract  N5ori-07850  with  the  Office  of  Naval 
Research. 

Two  reports  to  O.N.R.  on  this  study  will  treat  each  of  the  two  major 
phases  of  the  work. 

a.  "Theoretical  Background  of  Inertial  Navigation  for  Submarines",  which 
considers  the  physical  background  of  inertial  navigation  and  discusses 
idealized  systems  for  indication  of  position,  with  no  conclusions  drawn 
as  to  the  relative  suitability  of  the  different  systems  for  submarine 
operation. 

b.  "Characteristics  of  Systems  Feasible  for  Inertial  Navigation  of  Sub¬ 
marines",  which  considers  the  effect  of  imperfections  in  actual  instru¬ 
mentation,  e.g.,  "noise"  and  "drift",  on  system  performance,  leading 
to  recommendation  of  one  inertial  system  as  most  suitable  for  sub¬ 
marine  operation. 

Particular  mention  should  be  made  of  the  work  of  Mr.  Forrest  E.  Houston 
for  general  supervision  of  the  material  contained  in  the  report,  of  Mr.  John 
Hovorka  for  the  actual  writing  of  the  report,  and  the  Technical  Publications 
Division  of  Jackson  &'  Moreland  for  the  preparation  of  the  report. 


Cambridge,  Mass. 
March,  1951 


Walter  Wrigley 
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INTRODUCTION 


Inertial  guidance  of  submarines  can  be  considered  as  an  extension  of  nav¬ 
igation  techniques  to  self-contained  systems  that  operate  without  reference  to 
external  information,  save  for  initial  set-up  conditions.  It  is  sometimes  hu¬ 
morously  referred  to  as  "astronomy  in  a  closet".  In  this  respect  inertial  nav¬ 
igation  can  be  looked  on  as  a  precision  dead-reckoning  method  in  which  meas¬ 
urements  are  made  with  respect  to  inertial  space*  rather  than  with  respect  to 
the  water  surrounding  the  ship. 

Conventional  celestial  navigation  (see  any  text  on  navigation,  e.g., 
Dutton'1’**)  is  a  form  of  partial  inertial  navigation,  in  which  stars  represent 
the  inertial  space  references.  The  requirement  for  submerged  operation  pre¬ 
cludes  the  consideration  of  celestial  systems  in  this  report,  however.  In  full 
inertial  systems  gyroscopes  are  used  for  inertial  reference.  The  problem  of 
navigation  with  self-contained  systems  has  been  surveyed  by  Draper  and  his 
collaborators^),  and  the  basic  theoretical  background  of  inertial  guidance  has 
been  discussed  by  Wrigley''*'. 

The  first  part  of  this  report  treats  the  basic  physical  factors  that  are 
available  for  inertial  navigation.  This  work  is  based  primarily  on  that  of  ref¬ 
erence  (3),  but  particular  emphasis  is  placed  on  aspects  suitable  for  submarine 
operation.  A  discussion  of  the  Schuler  (84 -minute)  pendulum  tuning  condition, 
without  which  inertial  systems  would  scarcely  be  practical,  follows.  A  more 
complete  treatment  of  the  Schuler  tuning  condition  has  recently  appeared(29). 

The  final  part  of  the  report  treats  several  basic  systems  for  inertial  indication 
of  position.  This  treatment  is  idealized  in  that  ideal  components,  e.g.,  "drift¬ 
less"  gyro  units,  are  assumed,  for  ease  in  presenting  the  system  fundamentals. 
This  means  that  the  only  errors  discussed  in  this  report  are  those  due  to  pos¬ 
sible  initial  misalignment  of.  the  system.  Accordingly,  all  systems  treated  are 
considered  to  be  on  an  equal  basis,  and  no  attempt  is  made  in  this  report  to 
find  a  system  most  applicable  to  the  submarine  navigation  problem. 

In  the  next  and  final  report  these  potential  inertial  systems  will  be  exam¬ 
ined  from  the  practical  point  of  view,  i.e.,  from  the  point  of  view  of  the  effects 
of  actual  instrumentation  on  system  performance.  In  that  report  the  effects  of 
such  characteristics  as  integrator  drift  and  gyro  drift  will  be  considered,  and 
an  inertial  system  that  appears  most  suitable  for  submarine  navigation  will  be 
discussed  in  considerable  detail. 

It  will  be  noted  that  the  gyro  units  discussed  in  this  report  are  of  the  single - 
degree-of -freedom  type.  Although  it  is  fully  recognized  that  the  more  conven- 


*  Inertial  space  is  a  space  that  can  be  associated  with  the  "fixed  stars".  It  is  discussed  in 
detail  in  this  report. 

**  See  list  of  references  at  end  of  text. 


CONFIDENTIAL 


tional  two -degree-of -freedom  unit  could  just  as  well  be  used,  it  was  necessary 
to  restrict  the  discussion  to  one  type,  to  prevent  the  report  from  becoming  cum 
bersome,  and  because  the  former  type  of  unit  appears  to  be  more  widely  used 
in  present  systems  of  this  general  nature. 

The  text  contains  only  discussion  and  illustrations;  all  mathematical  work 
is  in  the  appended  derivations.  The  notation  used  in  the  derivations  is  of  the 
self-defining  type,  based  on  a  formulation  given  by  C.  S.  Draper,  Notes  on 
Instrument  Engineering,  M.I.T.  Instrumentation  Laboratory,  Cambridge, 

Mass.,  September,  1950. 
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NATURAL  PHENOMENA  AVAILABLE 
FOR  SUBMERGED  NAVIGATION 

Navigation  by  Radiation  Data 

Visual  navigation,  using  recognized  landmarks  and  associated  maps,  is 
a  special  class  of  navigation  using  radiation.  It  is  necessarily  subject  to 
artificial  and  natural  interference.  In  fact,  the  substitution  of  sound  or 
lower  frequency  electromagnetic  waves  for  light  waves  does  not  remove  the 
interference  limitation;  all  radiation  systems  have  this  in  common.  The 
interference-proof,  self-contained  navigation  systems  subsequently  described, 
therefore,  will  be  confined  to  the  observation  of  natural  phenomena  asso¬ 
ciated  with  the  Earth.  The  fundamental  kinematics  and  dynamics  of  such 
systems  are  the  subject  of  this  report. 

Utilization  of  Earth  Phenomena 

Navigation  by  Dead  Reckoning 

Three  quantities  available  for  submarine  navigation  by  dead  reckon¬ 
ing  are  the  following: 

1.  static  pressure, 

2.  dynamic  pressure,  and 

3.  temperature  of  the  surrounding  water. 

They  can  be  measured  directly  by  a  self-contained  navigation  system. 

From  these  data,  the  ship's  depth  and  its  velocity  (in  magnitude  and 
direction)  with  respect  to  the  water  can  be  determined.  However,  the 
direction  of  the  velocity  vector  with  respect  to  a  reference  direction 
fixed  in  the  Earth  cannot  be  determined  from  these  data.  Such  information 

*  requires,  among  other  things,  auxiliary  directional  indication,  e.g.,  a 
gyrocompass.  This  permits  the  indication  of  position  with  respect  to  the 
water,  although  not  with  respect  to  the  Earth.  Position  on  the  Earth  is 

•  established  when  local  ocean  current  data  are  used  to  modify  the  ship's 
velocity  information.  A  complete  self-contained  system  for  navigation  of 
a  submarine  anywhere  on  the  Earth  would  then  be  achieved. 

The  reliability  of  these  data  will  not  involve  merely  the  precision  with 
which  the  system  measures  them.  It  is  necessary  also  that  position  be 
uniquely  determined.  If  the  system  is  to  indicate  position  anywhere  on  the 
Earth,  there  must  be  a  one-to-one  correspondence  between  any  given  set 
of  data  and  a  position  on  the  Earth.  More  than  this,  Earth -fixed  data  which 
change  in  time—  in  particular,  ocean  current  data  —  must  be  known  at  all 
times.  This  problem  of  uniqueness  will,  in  general,  be  an  important  factor 
in  the  discussion  of  self-contained  navigation  systems. 


The  Earth's  Electric  Field  as  a  Basis  for  Navigation 

The  Earth  is  the  negative  surface  for  an  electrostatic  field  extending 
outward  into  space.  This  field  has  an  electrostatic  potential  gradient  of 
about  130  volts  per  meter  near  the  surface  of  the  Earth  and  diminishes 
with  altitude.  Very  little  further  information  applicable  to  navigation  is 
known.  Thus,  the  Earth's  electric  field  data  are  considered  inadequate 
as  the  basis  of  a  self-contained  navigation  system. 

The  Earth's  Magnetic  Field  as  a  Basis  for  Navigation 

The  Earth's  magnetic  field  has  been  used  in  the  guidance  of  vehicles 
for  centuries.  This  field  has  two  poles,  one  in  northern  Canada  and  one 
in  Antarctica,  and  effectively  acts  outside  the  mass  of  the  Earth  as  though 
it  were  due  to  either  a  spherical  magnet  or  an  imbedded  magnetic  dipole. (5)(6) 
The  field  is  normal  to  the  surface  of  the  Earth  at  the  magnetic  poles  and  tan¬ 
gent  to  the  Earth  at  the  magnetic  equator.  Accordingly,  the  magnetic  field 
does  not  have  a  unique  direction  at  any  given  point  on  the  Earth.  For  example, 
at  the  magnetic  equator  the  field  is  ideally  everywhere  parallel.  Thus  the 
magnetic  field  alone  will  not  uniquely  specify  position  on  the  Earth. 

When  auxiliary  equipment  is  used  to  establish  the  horizontal  plane,  the 
magnetic  field  offers  valuable  means  for  establishing  direction  in  the  plane. 
This  is  done  with  the  magnetic  compass.  Correction  for  the  angular  dif¬ 
ference  between  magnetic  north  and  geographic  north,  (declination  or  vari¬ 
ation)  utilizes  the  fact  that  the  properties  of  the  magnetic  field  are  quite  well 
known  for  all  of  the  surface  of  the  Earth,  except  in  the  polar  regions. 

The  total  intensity  of  the  magnetic  field  is  lowest  at  the  magnetic  equator 
and  increases  toward  the  magnetic  poles,  where  it  becomes  approximately 
twice  as  large  as  at  the  equator.  The  total  field  intensity  might  therefore 
serve  as  an  identification  for  position  on  the  Earth.  (7)  However,  such  a 
location  is  not  unique,  because  it  furnishes  only  a  locus  of  positions  having 
a  given  intensity.  (Such  a  locus  is  roughly  a  parallel  of  magnetic  latitude). 
Auxiliary  information  is  required  to  establish  a  unique  position  on  the  locus. 

When  an  electric  conductor  is  in  motion  relative  to  a  magnetic  field,  a 
potential  difference  proportional  to  this  velocity  is  induced  between  the  ends 
of  the  conductor.  (8)  it  has  often  been  suggested  that  the  current  associated 
with  such  a  potential  difference  could  be  used  to  indicate  ground  speed,  i.  e. , 
a  conductor  could  be  moved  so  as  to  cut  the  "lines  of  force"  of  the  Earth's 
magnetic  field.  The  fallacy  of  the  suggestion  lies  in  the  fact  that  the  electro¬ 
motive  force  required  to  drive  such  a  current  depends  on  the  net  rate  of  change 
of  magnetic  flux  through  a  closed  circuit.^)  In  any  circuit  that  can  be  carried 
in  a  vehicle,  the  Earth's  magnetic  flux  is  sufficiently  homogeneous  that  there 
is  no  net  change.*  Therefore  the  expected  current  is  always  zero,  and  no 
ground  speed  information  is  available  from  this  source. 


*  It  is  to  be  noted  that  in  a  closed  electric  circuit  moving  in  an  homogeneous  electromagnetic 
field,  exactly  as  much  (lux  enters  the  loop  as  leaves  it  in  a  given  time. 
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Magnetic  data  are  subject  to  interference  by  natural  causes,  notably 
magnetic  storms.  These  storms  appear  to  be  closely  associated  with  the 
aurora  and  sunspot  activity.  While  artificial  interference  with  magnetic 
data  cannot  be  ruled  out,  it  does  not  appear  likely  at  the  present  time  because 
of  the  large  amount  of  power  required  to  generate  any  effective  dynamic 
disturbance.  A  static  disturbance  would  be  encountered  in  the  shielding 
effect  of  the  hull  of  the  ship  carrying  the  navigation  system.  Accordingly, 
in  a  submerged  submarine,  the  steel  pressure  hull  would  effectively  shield 
out  magnetic  field  data. 

The  Earth’s  Gravity  Field  as  a  Basis  for  Guidance 

According  to  Newton's  law  of  gravitation,  every  mass  particle  in  the 
universe  attracts  every  other  particle  with  a  force  proportional  to  the  product 
of  their  masses  divided  by  the  square  of  the  distance  between  them.  Any 
massive  body,  such  as  the  Earth,  can  be  considered  as  having  associated 
with  it  a  gravitational  field.  Such  a  field  is  essentially  radial  to  the  center 
of  the  Earth  and  has  no  poles  on  the  Earth's  surface.  The  possibility  there¬ 
fore  exists  for  the  direction  of  this  field  at  any  given  point  on  the  Earth  to 
serve  as  a  unique  identification  of  the  position  of  that  point.  The  direction 
of  the  gravitational  field  could  be  readily  determined  by  a  plumb  bob,  if  the 
Earth  were  non-rotating.  The  daily  rotation  of  the  Earth,  however,  produces 
a  centrifugal  acceleration  that  also  affects  the  plumb  bob.  This  is  due  to  a 
basic  law  of  physics,  namely  the  principle  of  equivalence  in  the  general 
theory  of  relativity,  which  states  that  gravitational  mass  is  equivalent  to 
inertial  mass.  Because  of  this,  it  is  impossible  to  distinguish  directly 
between  gravitational  effects  and  linear  accelerations/*®)  The  vector 
resultant  of  the  Earth's  attracting  gravitational  field  and  the  centrifugal 
force  field  due  to  the  Earth's  daily  rotation  is  defined  as  the  Earth's  gravity 
field  (see  Fig.  1).  The  direction  of  this  essentially  radial  gravity  field  is 
unique  at  any  given  point  on  the  Earth.  The  secular  (time)  variation  in  the 
direction  of  the  gravity  fiqld  at  a  given  point  is  caused  mainly  by  tidal  effects, 
and  is  less  than  0.05  microradian.  (**)  For  this  reason  the  direction  of  the 
gravity  f-ield  is  a  reliable,  unique  characteristic  of  any  given  point  on  the 
Earth.  Furthermore,  it  is  essentially  impossible  to  interfere  with  gravity 
effects.  Accordingly,  the  Earth's  gravity  field  appears  to  be  the  most 
promising  source  of  data  for  a  self-contained  navigation  system. 

The  Earth  itself  is  also  subject  to  the  equivalence  of  gravitational  and 
inertial  effects.  If  the  Earth  were  a  non-rotating  fluid  mass,  it  would  form 
a  sphere  under  the  influence  of  its  own  gravitational  attraction.  Because  of 
the  Earth's  structure  and  daily  rotation,  however,  the  associated  centrifugal 
force  field  causes  the  Earth  to  bulge  at  the  equator,  giving  rise  to  a  nearly 
ellipsoidal  or  spheroidal  shape.  The  figure  that  a  fluid  body  with  the  mass 
distribution  and  daily  rotation  of  the  Earth  would  have  is  defined  as  the 
geoid,  a  "square  shouldered"  ellipsoid  with  a  slightly  undulated  surface. 

The  surface  of  the  geoid  is  represented  by  mean  sea  level.  Variations  in 
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the  elevation  of  the  geoid  relative  to  the  closest  reference  ellipsoid  are 
approximately  one  percent  of  the  topographic  variations  in  elevation. 

Since  the  geoid  does  not  form  an  analytical  surface,  attempts  are  made 
to  represent  the  geoid  by  an  associated  reference  ellipsoid  for  the  basis  of 
geodetic  measurem'ents  and  map  data.  The  Hayford  Spheroid  of  1909,  which 
is  the  international  reference  ellipsoid,  *  has  an  equatorial  radius  of 
6,378,388  meters  and  a  polar  semi-diameter  of  6,356,912  meters,  which 
gives  an  ellipticity**  of  1/297.  (11)(12)(13) 


Figure  1.  Gravitational  Field  Intensity  and  Gravity  Field  Intensity  Ithe  Vertical). 

The  surface  of  the  geoid  is  everywhere  an  equipotential  surface  of  the 
Earth's  gravity  field.  The  direction  of  the  gradient  of  the  gravity  potential 
at  the  surface  of  the  geoid  is  defined  as  the  vertical.  It  is  the  direction  of 
the  specific  force  of  gravity***  as  indicated  by  a  plumb  bob  with  its  base 
stationary  with  respect  to  the  Earth  (Fig.  1).  The  specific  force  of  gravity 
(denoted  by  the  vector  g)  is  the  vector  resultant  of  the  gravitational  specific 
force  (denoted  by  the  vector  G)  and  the  centrifugal  specific  force  associated 
with  daily  rotation. 

Because  the  geoid  does  not  have  a  consistently  smooth  surface,  the 
vertical  is  generally  not  parallel  to  the  normal  to  the  reference  ellipsoid 

*  The  Clarke  Snherold  of  1866  Is  the  basis  for  North  American  triangulations. 

**  Ellipticity  is  defined  as  the  ratio  of  the  difference  between  equatorial  radius  and  the  polar 
semi-diameter  to  the  equatorial  radius. 

***  The  specific  force  of  gravity  is  defined  as  the  force  of  gravity  acting  on  a  unit  mass.  It  is 
identical  with  the  gravity  field  Intensity.  Gravitational  specific  force  is  identical  with  the  gravi¬ 
tational  field  intensity.  Centrifugal  specific  force,  which  is  a  reaction  force,  is  equal  in  magni¬ 
tude  but  opposite  in  direction  to  the  associated  centripetal  acceleration.  Specific  force,  expressed 
in  the  dimensions  of  force  per  unit  mass,  is  numerically  and  dimensionally  equal  to  acceleration. 
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at  the  same  position.  The  angle  between  the  vertical  and  the  ellipsoid 
normal  is  called  the  deflection  of  the  vertical.  It  is  generally  less  than 
0.30  milliradian.  In  general,  the  largest  marine  deflections  occur  in  the 
regions  of  rapid  changes  in  ocean  depth,  such  as  near  deeps  or  islands 
that  rise  abruptly  from  the  ocean  bottom. 

The  position  of  a  point  on  the  Earth's  surface  can  be  determined 
either  by  geodetic  measurements,  on  which  maps  are  based,  or  by  astro¬ 
nomical  position  measurements,  which  utilize  the  vertical.  (The  measure¬ 
ment  of  astronomical  position  is  subsequently  discussed  in  this  report.) 
Geodetic  and  astronomical  position  data  will  not  agree,  in  general,  because 
of  the  above-mentioned  deflection  of  the  vertical.  However,  position  data 
determined  by  vertical  indication  and  corresponding  map  data  should 
agree  in  the  case  of  small  land  masses  such  as  islands  whose  position 
have  been  mapped  from  astronomical  data. 

In  general,  variations  in  the  deflection  of  the  vertical  can  be  con¬ 
sidered  to  be  quasi-static  at  submarine  speeds,  and  can  be  expected  to 
have  no  dynamic  effect  on  the  operation  of  submarine  inertial  navigation 
systems. 
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GRAVITY  FIELD  NAVIGATION 
Astronomical  Position 

The  uniqueness  of  the  vertical  at  any  given  point  on  the  earth  is  the  basis 
for  the  astronomical  position(H)  of  the  point.  Astromonical  position  is  an 
angle  centered  in  the  Earth  and  has  two  components,  namely, 

1.  Astronomical  latitude,  the  angle  between  the  equatorial  plane  and 
the  vertical  at  the  located  point.  This  quantity  is  unique  because 
the  Earth's  polar  axis  has  a  unique  spatial  direction. 

2.  Astronomical  longitude,  the  angle  between  the  plane  through  the 
Earth's  polar  axis  containing  an  arbitrary  vertical  (e.g. ,  that 
at  Greenwich)  and  the  plane  through  the  polar  axis  containing  the 
located  vertical.  This  quantity  itself  is  not  unique  because  of  its 
arbitrary  reference.  However,  difference  of  longitude  is  unique. 

Inertial  Navigation:  Definition 

Inertial  navigation  of  a  vehicle  is  defined  as  a  method  of  locating  the 
vehicle's  position  with  respect  to  inertial  space,  using  gravity  field  and 
acceleration  data  and  Newton's  laws~of  motion.  *  The  concepts  stated  and 
implied  in  this  definition  will  now  be  discussed. 

Newton's  first  law  of  motion  states  that  a  body  suffering  no  external 
forces  is  not  accelerated,  and  therefore  moves  with  constant  velocity,  i.e. . 
constant  speed  in  a  straight  line.  The  coordinate  system  to  which  this 
velocity  is  referred  is  called  an  inertial  space.  (10)  (14)  (15)  Since,  in 
particular,  a  centripetal  force  is  not  to  act  on  the  body,  the  coordinate 
system  must  be  non-rotating.  It  can  only  move  (with  respect  to  some  other 
coordinate  system)  along  a  straight  line  with  constant  velocity.  Strictly 
speaking,  it  is  impossible  to  distinguish  in  an  absolute  sense  between  two 
such  coordinate  systems,  although  Newtonian  mechanics  associates  an 
absolute  inertial  space  with  the  average  positions  of  the  so-called  "fixed 
stars."  However,  special  relativity  theory  associates  the  idea  of  an  inertial 
space  with  any  coordinate  system  moving  at  constant  velocity  relative  to  the 
material  bodies  of  the  universe.  Any  space  fulfilling  this  requirement  can  be 
chosen  as  a  reference  inertial  space.  Because  of  the  principle  of  equivalence 
in  the  general  theory  of  relativity,  it  is  necessary  to  include  gravitational 
effects  in  determining  an  inertial  space. (10)(15) 

Equivalence  of  Gravitational  and  Inertial  Effects 

One  of  the  largest  fundamental  problems  associated  with  inertial  navi¬ 
gation  is  the  indication  of  the  vertical  (e.  g. ,  the  Earth's  gravity  field  effect 
on  a  plumb  bob;  a  gravitational**  effect)  inside  a  self-contained  vehicle  that 


See  any  standard  physical  mechanics  text  for  reference. 

For  this  illustration  the  difference  between  gravity  and  gravitation  is  ignored. 
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is  suffering  additional  linear  accelerations  (inertial  effects).  The  principle  of 
equivalence  in  the  theory  of  general  relativity  states  that  gravitational  and  in¬ 
ertial  effects  are  intrinsically  indistinguishable.  A  simple  example  (16)  Can 
illustrate  this  point:  a  man  stands  in  an  elevator,  at  rest  with  respect  to  the 
Earth,  his  feet  on  a  platform  scale.  The  scale  reads  his  weight,  i.e.,  the  Earth 

pulls  on  him  with  a  force  equal  to  his 
weight;  but  he  does  not  accelerate,  be¬ 
cause  the  scale  pushes  up  on  his  feet 
with  a  force  of  the  same  magnitude  as 
his  weight.  Now  suppose  the  elevator 
accelerates  upward,  propelled  by  an 
external  force.  The  scale  now  reads 
the  inertia  reaction  force  associated 
with  this  acceleration  in  addition  to 
his  weight.  He  observes  that  he 
weighs  more;  that  is,  it  is  as  if  the 
Earth's  gravitational  field  had  in¬ 
creased.  If  the  external  force  is 
now  removed  and  the  elevator  falls 
freely  in  the  Earth's  gravitational 
field,  man,  scale  and  elevator  have 
the  same  downward  acceleration. 

The  scale  reads  zero,  and  it  is  as  if 
the  Earth's  gravitational  field  were 
zero,  as  far  as  the  man  in  the  ele¬ 
vator  is  concerned.  Thus,  from  with¬ 
in  the  elevator,  an  observer  cannot 
distinguish  gravitational  from  acceler¬ 
ation  effects.  -V- 


The  X,  Y,  Z  axes  in  tiicsfi  system?  are  mutually  orthogonal 

Inertial  Space  Reference  Axes  -  Centered  at  center  of  Earth 

Zj  Parallel  to  polar  axis 

Xj.  Y«  Arbitrarily  located  in  equatorial  plane 

Non-rotating  with  respect  to  "fixed  stars” 

Earth  Reference  Axes  -  Centered  at  center  of  Earth 

!£g  Parallel  to  polar  axis 

YE  Arbitrarily  located  In  equatorial  plane 
Non-rotating  with  respect  to  the  Earth 

Earth  reference  axes  make  one  rotation  about  the  polar  axis  with 
respect  to  inertial  space  reference  axes  In  one  sidereal  day. 


Figure  2.  Geometrical  Relationships  between 
Earth  Reference  Axes  and  Inertial  Space 
Reference  Axes. 

Coordinate  Systems 

For  navigation  purposes  it  would  be  desirable  if  acceleration  could  be 
measured  directly  with  respect  to  the  Earth.  Unfortunately  for  the  designer, 
the  acceleration-sensitive  equipment  to  be  discussed  makes  its  measurements 
with  respect  to  inertial  space,  for  which  reason  it  is  necessary  to  evaluate 
the  several  acceleration  components  that  make  up  these  measurements.  This 
can  be  facilitated  by  a  proper  choice  of  coordinate  systems,  of  which  four 
adequately  serve  the  needs  of  basic  navigation: 


1.  Inertial  space  reference  system 

2.  Earth  reference  system 

3.  Position  reference  system 

4.  Position  system 


Each  of  these  is  subsequently  discussed  in  detail. 
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Inertial  Space  Reference  System 

The  inertial  space  reierence  system  must  be  one  in  which  Newton's  law 
of  inertia  is  valid.  This  system  can  be  most  conveniently  centered  at  the 
center  of  the  Earth,  which  is  allowable  because  the  Earth  is  in  an  orbital,  or 
free-fall,  condition  about  the  Sun*.  The  axes  of  this  inertial  space  reference 
system  must,  be  non-rotating  with  respect  to  the  "fixed  stars".  The  Earth- 
centered  inertial  space  reference  system ’is  illustrated  in  fig.  2.  Here,  one 
of  the  three  orthogonal  axes,  marked  Zp  is  along  the  Earth's  polar  axis,  and 

the  other  two,  X^  and  Yp  are  arbitrarily  located  in  the  plane  of  the  Earth's 

equator,  and  are  non-rotating  with  respect  to  the  "fixed  stars". 

Earth  Reference  System 

Earth  reference  axes  are  also  most  conveniently  centered  at  the  center  of 
t?  Earth  and  are  non-rotating  with  respect  to  the  Earth.  As  shown  in  Fig.  2, 
o..e  reference  axis,  Z^,  is  selected  to  be  coincident  with  the  Earth's  polar 

axis,  and  the  other  two  axes,  X„  and  Y_,  are  arbitrarily  located  in  the  plane 

tj  hi 

of  the  equator.  Earth  reference  axes  make  one  rotation  about  the  polar  axis 
with  respect  to  the  inertial  space  reference  axes  in  one  sidereal  day. 


*  In  an  orbit,  the  central  gravitational  field  exactly  cancels  the  centrifugal  field,  thus  leading  to 
zero  net  acceleration. 


INDICATED  POSITION  TRUE  POSITION 


■  tUAC*****  *«ft’  <V’I  <• 

MU'Mwm'iOA  KMtvn 

*»  ren 


CORRECTION  TO  INDICATED 
VERTICAL  WITH  RESPECT  TO 
TRUE  VERTICAL  (A  VECTOR 
ANGLE) 

CC)(v„.,(tl)3  KC)V|(U) 


Itio  X,  v  Z  axes  in  these  systems  are  mutually  orthogonal. 

True  Position  Axes  -  Centered  at  Center  of  Earth. 

Zp|  Parallel  to  true  vertical. 

Xpt  and  Ypi  are  parallel  to  the  direction  ot  the  v/ue 

north  and  normal  to  the  meridian  plane  through  the 
true  position  respectively. 

Indicated  Position  Axes  -  Centered  at  Center  of  Earth. 

ZPl  Para,lcl  t0  indicated  vertical  (l  e.,  vertical  at 
Indicated  position.) 

Ypt  and  Xpt  are  rotated  from  Yp<  and  Xp<  respect¬ 
ively  through  the  vector  angle  ((C)V)(t  ^  plus  any 
difference  between  a  true  and  Indicated  astmuth. 

A  vector  angle  is  characterized  by  the  direction  o(  the  axis 
about whichttrepres'ntsa  rotation  aswellasby  a  magnitude. 


Figure  3.  True  Position  and  Indicated  Position  Coordinate  Systems. 
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Position  Reference  System 

The  position  reference  coordinate  system  for  navigation  purposes  will  be 
taken  as  an  Earth-centered  system,  the  Z-axis  coincident  with  the  true  verti¬ 
cal  at  the  vehicle's  present  position.  The  Y-axis  is  arbitrarily  taken  as  the 
normal  to  the  meridian  plane  through  the  true  position,  while  the  X-axis  is 
parallel  to  the  direction  of  true  north.  These  axes  are  shown  as  X  .,  Y  , 
and  Zpt  in  Fig.  3. 

Position  System 

The  position  coordinate  system  for  navigation  purposes  will  also  be  taken 
as  an  Earth-centered  system,  the  Z-axis  in  this  case  being  coincident  with 
the  indicated  vertical  at  the  vehicle's  present  position.  Indicated  axes  (Xp., 

Yp.,  Zp.,  in  Fig.  3)  are  identical  with  true  axes  when  the  indicated  position 

and  indicated  azimuth  reference*  of  the  indicated  position  system  are  iden¬ 
tical  with  the  true  position  and  true  azimuth  reference.  When  the  two  posi¬ 
tions  and  horizontal  references  are  not  the  same,  the  indicated  axes  are 
rotated  from  the  true  axes  by  the  vector  sum  of  two  angles(^): 

1.  The  angle  between  the  true  vertical  and  the  indicated  vertical 
([(C)V](M)  in  Fig.  3) 

2.  The  angle  between  the  true  horizontal  reference  and  the  indicated 
horizontal  reference. 

Practical  Coordinate  Systems 

The  coordinate  systems  illustrated  in  Figs.  2  and  3  serve  well  for  the 
theoretical  background  because  they  are  geocentric.  In  any  actual  navigation 
system,  however,  the  coordinate  origin  will  not  be  at  the  Earth's  center, 
but  in  the  vehicle,  and  the  coordinate  axes  will  be  parallel  to  the  earth 
reference  axes  or  inertial  axes.  Moreover,  additional  coordinate  systems 
will  be  required,  some  fixed  in  the  vehicle,  some  in  the  indicating  equip¬ 
ment,  etc.,  to  account  for  the  actual  conditions  encountered  in  the  solution 
of  a  navigation  problem.  Such  conditions  complicate  the  problems  of  the 
system  designer,  but  they  do  not  invalidate  the  conclusions  reached  in  this 
paper  for  the  idealized  theoretical  situation. 

Theory  of  Vertical  indicators 

Specific  Forces  Associated  with  Gravity  Field  Navigation 

The  simplest  vertical  indicator  is  a  plumb-bob  or  pendulum.  In  general, 
a  pendulous  element  is  defined  as  a  system  free  to  rotate  about  a  point,  so 
that  it  assumes  a  preferred  orientation  under  the  influence  of  impressed 
external  forces.  If  the  suspension-point  is  not  accelerated  with  respect  to 
the  Earth,  a  pendulous  element  can  provide  a  precise  indication  of  the  true 


*  Indicated  azimuth  represents  a  rotation  about  the  indicated  vertical. 
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vertical,  i.e.,  the  direction  of  gravity.  If  the  suspension-point  is  accelerated 
with  respect  to  the  Earth,  the  pendulous  element  tends  to  indicate  the  apparent 
vertical,  i.e.,  the  direction  of  the  resultant  specific  force. 


When  the  pivot  of  a  pendulous  element  is  accelerated,  a  reaction  force  acts 
on  the  center  of  mass  of  the  pendulous  element.  This  force  produces  a  torque 
that  causes  the  pendulous  element  to  "lag  behind"  the  pivot.  To  study  the  effects 
of  an  acceleration  without  introducing  the  effect  of  the  mass  of  the  pendulous 
element,  it  is  convenient  to  discuss  the  motion  in  terms  of  force  per  unit  mass, 

i.e.,  in  terms  of  the  inertia  reaction 
specific  force  associated  with  the 
acceleration.  An  inertia  reaction 
specific  force  has  the  magnitude  of 
its  associated  acceleration  (although 
it  will  be  expressed  in  force-per- 
mass  units  rather  than  in  accelera¬ 
tion  units)  and  is  opposite  in  direc¬ 
tion.  Its  direction  is  the  one  that  a 
pendulous  element  seeks  (18);  and 
its  projection  on  the  input  axis  of  a 
linear  accelerometer  is  the  direc¬ 
tion  in  which  the  accelerometer 
element  is  urged.  As  mentioned 
earlier,  the  specific  gravitational 
force  is  identical  with  the  conven¬ 
tional  gravitational  field  intensity. 
Figure  4  is  a  force  diagram  for  an 
accelerated  pendulous  element  in  a 
gravitational  field. 


Figure  4. 


Specific  Forces  Acting  on  a  Pendulous  Element.  To  obtain  the  Specific  forces 

acting  on  a  pendulous  element  during 

a  navigation  problem  it  is  necessary  to  know  the  associated  accelerations. 

Page,  in  a  reference(^),  shows  that  the  linear  acceleration  of  a  point  with  res¬ 
pect  to  a  reference  coordinate  system  may  comprise  five  components*  with 
respect  to  a  coordinate  system  that  is  moving  relative  to  the  reference  coordin¬ 
ate  system**.  These  components  are  listed  in  column  A  at  the  left  side  of  Table  1. 
The  table  lists  also  the  results  of  Derivation  1  (appended)  in  columns  B  and  C, 
giving  the  corresponding  acceleration  components  of: 


1.  A  vehicle  whose  true  position  coordinates  move  with  respect  to  an 
Earth  reference  coordinate  system;  acceleration  components  are  taken 
with  resnect  to  the  Earth. 


*  Wrigley  has  elsewhere  expanded  this  list  to  include  simultaneous  relationships  between 
several  coordinate  systems. 

**  It  Is  to  be  noted  that  the  reference  coordinate  system  may  also  be  moving  with  respect  to  in¬ 
ertial  space,  but  such  motion  is  non-essential  to  this  discussion. 


ACCELERATION  COMPONENTS  ASSOCIATED  WITH  SPECIFIC  FORCES 
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2.  A  vehicle  in  a  true  position  system  stationary  with  respect  to  the  Earth 
reference  coordinate  system;  acceleration  components  are  taken  with 
respect  to  an  inertial  coordinate  system. 

As  far  as  marine  navigation  is  concerned,  only  two  components  in  columns 
B  and  C  of  Table  1  are  significant: 

1.  The  tangential  (horizontal)  acceleration  associated  with  the  angular 
acceleration  of  the  true  vertical  as  the  vehicle  accelerates  over  the 
Earth's  surface. 

2.  The  centripetal  acceleration  parallel  to  the  projection  of  the 
Earth-radius  on  the  equatorial  plane,  due  to  the  Earth's  daily 
rotation;  this,  added  to  the  gravitational  field  intensity,  gives 
the  gravity  field  intensity. 

Geocentric  Rotations 

The  discussion  thus  far  has  made  repeated  use  of  the  concept  of  a  rotating 
true  vertical  associated  with  the  present  position  of  the  vehicle  on  the  Earth's 
surface.  This  rotation  of  the  vertical  is  geocentric,*  and  gives  rise  to  two  im¬ 
portant  acceleration  components  of  the  vehicle,  as  stated  above.  The  concept 
is  fundamental;  essentially,  all  gravity  field  navigation  systems  are  vertical 
indicators.  The  direct  applicability  of  vertical  indication  to  navigation  stems 
from  the  already-mentioned  fact  that  the  vehicle's  astronomical  position  can 
be  expressed  in  terms  of  verticals  taken  at  various  points  on  the  Earth.  The 
application  of  a  vertical  indicator  to  navigation  admits  of  what  appear  super¬ 
ficially  to  be  two  distinct  points  of  view.  In  one  case,  the  comparison  of  the 
indicated  vertical  with  a  reference  vertical  is  used  to  get  the  desired  latitude 
or  longitude.  The  indicated  vertical  is  then  a  datum  involved  in  determining 
present  position.  According  to  the  second  point  of  view,  the  geocentric  angu¬ 
lar  displacement  of  the  vehicle  is  derived  from  double  integration  of  the  tan¬ 
gential  acceleration  measured  by  the  navigation  apparatus,  and  the  vertical  in¬ 
dicator  plays  the  subsidiary  but  entirely  necessary  role  of  orienting  the  appa¬ 
ratus  so  that  it  will  measure  tangential  accelerations  only,  i.e.,  it  effectively 
provides  a  stable  platform  to  eliminate  the  effect  of  gravity  in  the  acceleration¬ 
measuring  equipment.  In  either  case,  the  basic  inseparability  of  tangential  and 
vertical  accelerations  requires  the  presence  of  the  vertical  indicator  as  such. 
These  ideas  will  be  discussed  at  length  later. 

Fundamentals  of  the  Schuler  (84-Minute)  Pendulum 

When  a  vehicle  moves  over  the  nearly  spherical  earth,  the  vertical 
associated  with  its  instantaneous  position  rotates  approximately  geo¬ 
centrically  with  respect  to  the  Earth.  On  a  ship  moving  at  constant 
sneed  along  a  great  circle  over  the  Earth,  a  pendulum  would  have  to 
rotate  with  respect  to  the  Earth  at  an  effectively  constant  angular  velo¬ 
city  in  order  to  indicate  the  vertical  continuously.  Once  the  transient 

*  Strictly  sneaking,  this  rotation  is  not  truly  geocentric,  since  the  Earth  is  not  actually  SDherical. 
However,  in  the  subsequent  text,  this  rotation  will  be  referred  to  as  "essentially  geocentric". 
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stages  were  over,  no  torque  would  be  required  to  keep  the  pendulum  along 
the  vertical  during  this  rotation.  On  the  other  hand,  if  the  ship  were  to  ac¬ 
celerate,  by  changing  either  speed  or  heading,  a  torque  would  be  required 
to  maintain  the  pendulum  on  the  vertical,  which  in  turn  is  geocentrically 
accelerating. 

Schuler  pointed  out  in  1923  that  accurate  indication  of  the  vertical 
during  periods  of  acceleration  could  be  realized  by  suitable  tuning  of  the 
natural  period  of  the  pendulum.  Whenever  the  pivot  of  a  pendulum  is  accele¬ 
rated,  the  center  of  mass  of  the  pendulum  tends  to  "lag  behind”  the  pivot 
with  respect  to  inertial  space.  At  the  same  time,  this  acceleration  causes 
the  true  vertical,  associated  with  the  pivot,  to  accelerate  geocentrically  with 
respect  to  inertial  space.  These  considerations  lead  to  the  following  obser¬ 
vation: 


If  a  pendulum  initially  hangs  vertically,  it  will  remain 
along  the  vertical  if  its  angular  acceleration  about  its 
pivot  equals  the  geocentric  angular  acceleration  of  the 
vertical.  (The  accelerations  are  with  respect  to  inertial 
space,  as  stated.) 

These  two  accelerations  become  equal  for  a  distributed-mass  pendulum 
when  the  ratio  of  the  pivot-center  of  mass  separation  to  the  square  of  the 
radius  of  gyration  of  the  pendulum  equals  the  reciprocal  of  the  radius  of  the 
Earth.  When  operation  is  in  the  Earth's  gravity  field,  and  the  pendulum  is 
undamped,  this  condition  gives  a  pendulum  with  a  natural  period  of  approxi¬ 
mately  eighty-four  minutes;  hence  the  term  84-minute  pendulum.  For  a 
shorter-period  pendulum  the  rotation-producing  torque  on  the  pendulous  • 
element  about  its  pivot  is  stronger,  so  the  pendulum  lags  the  true  vertical. 
For  a  longer-period  pendulum  the  torque  is  weaker  and  the  pendulum  leads 
the  true  vertical. 

Derivation  2  presents  a  short  mathematical  treatment  of  the  motion  of  a 
pendulum  under  the  influence  of  a  specific  force.  The  dependent  variable  in 
the  derivation  is  the  correction  that  must  be  applied  to  the  pendulum  in  order 
to  align  the  indicated  vertical,  which  is  established  by  the  pendulum,  with  the 
true  vertical  (see  Fig.  2-1  for  definitions  and  illustration  of  the  quantities 
involved).  This  variable  permits  direct  expression  of  the  performance  of  a 
pendulum  as  an  indicator  of  the  vertical. 

It  is  to  be  noted,  from  Fig.  2-1,  that  the  horizontal  inertia  reaction 
specific  force  is  equal  to  the  negative  of  the  product  of  the  radius  of  the 
Earth  and  the  angular  acceleration  of  the  true  vertical.  In  other  words,  as¬ 
suming  that  motion  over  the  Earth  involves  operation  over  an  essentially 
spherical  surface,  any  linear  acceleration  of  a  point  over  the  surface  is  in¬ 
evitably  accompanied  by  an  essentially  geocentric  angular  acceleration  of  the 
vertical  representing  that  point.  Because  of  the  aforementioned  principle  of 
equivalence  of  gravitational  and  inertial  masses,  and  assuming  operation  over 
a  spherical  surface,  the  orientation  of  the  vertical  (gravitational  effect)  and 
its  angular  acceleration  (inertial  effect)  are  inextricably  bound  together. 
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Clemens  and  his  collaborators  have  made  an  extensive  study  of  the 
characteristics  of  the  Earth-radius  pendulum,  i.e.,  the  Schuler-tuned  simple 
(concentrated  mass)  pendulum. 

An  84-minute  period  is  very  long  compared  to  the  periods  of  physical 
pendulums  ordinarily  encountered.  In  fact,  it  is  improbable  that  a  simple  * 

(concentrated  mass)  pendulum  or  a  distributed-mass  pendulum  could  be  con¬ 
structed  with  this  period, (20)  although  the  84-minute  gyropendulum  just  fails 
of  being  realizable.  Howe.er,  in  the  practical  case,  a  condition  analogous  to 
the  Schuler  tuning  of  a  pendulum  can  be  realized  in  a  particular  kind  of  closed- 
loop  system  in  which  a  short-period  pendulum  tracks  the  specific  forces.  A 
platform,  controlled  on  the  basis  of  the  pendulum  output  data,  then  becomes 
an  equivalent  Schuler  pendulum.  This  will  be  discussed  at  length  later. 

It  is  to  be  noted  that  there  is  a  definite  functional  relationship  between  the 
torque  applied  to  a  pendulum  about  horizontal  axes  and  the  linear  acceleration 
of  the  unit.  In  a  practical  vertical-indicating  system  this  torque-acceleration 
relationship  may  not  hold  due  to  added  properties  of  the  system,  as  is  demon¬ 
strated  in  the  subsequent  discussion.  Such  a  system  is  then  not  an  actual 
Schuler  pendulum  but  merely  exhibits  the  Schuler  tuning  condition  as  one  of  its 
characteristics. 

An  implicit  corollary  of  Schuler  pendulum  theory  is  that  the  direction  with 
respect  to  inertial  space  of  the  resultant  specific  force  remains  fixed  when  the 
pivot  of  the  pendulum  undergoes  a  harmonic  oscillation  with  an  84 -minute 
period.  In  general,  this  means  that  the  specific -force -tracking  element  in  a 
Schuler-tuned  inertial  navigation  system  will  be  insensitive  to  an  oscillatory 
disturbance  having  an  eightv-four  minute  period. 

The  Pendulum  as  a  Vertical  Indicator:  Summary 

The  following  ideas  are  deduced  from  the  foregoing: 

1.  A  pendulum  of  arbitrary  period  tracks,  i.e. ,  tends  to  align  itself  with 
the  apparent  vertical.  If  the  suspension-point  of  the  pendulum  is  ac¬ 
celerated,  the  true  vertical  is  then  not  available  as  a.  datum. 

2.  A  pendulum  with  Schuler  tuning  effectively  tracks  the  true  vertical  at 
all  times. 

3.  The  Schuler  pendulum  is  per  se  physically  unrealizable,  chiefiy  due  ,  I 

to  a  prohibitively  small  pivot-center  of  mass  separation. 

i 

4.  A  vertical  indicating  system  is  conceivable  which,  while  not  a  pendu¬ 
lum,  utilizes  the  Schuler  tuning  idea,  so  that  a  torque  caused  by  the  j 

shin's  acceleration  acting  on  the  system  will  be  absorbed  by  the  in¬ 
ertia  reaction  torque  required  to  keep  a  reference  line  in  the  system 

on  the  vertical;  this  is  the  fundamental  objective  of  Schuler  tuning. 


I 
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Smoothing,  Solution  Time  and  Forced  Dynamic  Error 

In  general,  no  data  received  by  an  actual  system  are  perfect,  and  the 
equipment  will  generate  some  uncertainty  in  its  output.  For  this  reason  it  is 
necessary  to  incorporate  some  degree  of  smoothing  into  any  navigation  sys¬ 
tem  if  the  system  is  to  be  operated  over  extended  periods  of  time.  Smoothing, 
which  has  been  extensively  investigated  in  the  fields  of  communications,  auto¬ 
matic  control,  radar,  and  fire-control  (22) (23) (24)(25)  involves  the  suppression 
of  the  roughness  components  of  an  input.  Since  roughness,  or  noise,  generally 
appears  among  high-frequency  Fourier  components  of  a  signal,  smoothing 
means  attenuation  of  the  high-frequency  end  of  the  input  spectrum,  including 
legitimate  high-frequency  components  of  the  input  signal. 

This  high-frequency  suppression  "rounds  off"  the  response  of  a  system 
to  sudden  changes  in  the  input,  resulting  in  a  delay  in  the  system's  achieving 
the  state  associated  with  new  input  conditions.  One  convenient  measure  of 
such  a  delay  is  solution  time,  defined  for  this  report  as  the  time  for  the  sys¬ 
tem  to  achieve  ninety -five  percent  of  the  change  associated  with  a  given  input 
change  (24)*.  Solution  time  is  important  in  determining  the  time  required  to 
remove  transient  effects  associated  with  either  the  initial  phase  of  a  problem 
or  sudden  changes  in  the  input  quantity,  i.e.,  the  time  required  for  the  system 
to  generate  a  "solution".  In  general,  greater  smoothing  means  a  longer  solu¬ 
tion  time. 

When  the  input  changes  sufficiently  rapidly  with  respect  to  system  para¬ 
meters,  over  an  extended  period  of  time,  so  that  the  output -input  relationship 
is  a  function  of  the  change,  the  delaying  action  of  smoothing  leads  to  a  forced 
dynamic  error.  Increased  smoothing  increases  the  forced  dynamic  error. 

The  major  problem  facing  the  designer  of  an  operating  system  is  to 
determine  the  best  compromise  between  the  smoothing  required  to  furnish 
an  acceptable  answer  to  the  engineering  problem,  on  the  one  hand,  and  on 
the  other,  the  associated  delay  and  forced  dynamic  error  allowable  in 
achieving  the  solution. 

Pendulums  and  Accelerometers 

Basically,  pendulums  and  accelerometers  both  derive  data  from  the 
resultant  specific  force  acting  on  the  system.  The  essential  geometrical 
features  associated  with  operation  of  a  pendulum  are  illustrated  in  Fig.  5. 

A  two-degree-of -freedom  pendulum  tracks  the  resultant  specific  force  in 
such  a  manner  as  to  try  to  align  a  tracking  line,  i.e. ,  its  pivot-center  of 
mass  separation  vector,  with  the  direction  of  the  resultant  specific  force. 

The  outputs  of  such  a  unit  in  a  vertical  indicating  system  are  the  angle 


*  This  is  equivalent  to  the  time  required  to  realize  ninety-five  percent  of  the  area  under  the 
smoothing  function^). 


Single-Degree-of-Freedom  Pendulous  Elements  with  Mutually  Perpendicular  Output 
Axes  to  One  Ideal  Two-Degree-of-Freedom  Pendulous  Element. 
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* 


between  the  indicated  vertical  and  the  pendulum  tracking  line,  and  an  orien¬ 
tation*  of  the  tracking  line  as  a  rotation  about  the  indicated  vertical  with  re¬ 
spect  to  some  reference  in  the  indicated  horizontal  plane.  A  single -degree - 


of-freedom  pendulum  tracks  the  resultant  specific  force  in  such  a  manner  as 
to  try  to  align  its  tracking  line  with  the  direction  of  the  projection  of  the  re- 
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Figure  6.  Specific  Force  Vector  Diagram  for 
Two  Ideal  Single-Degree-of-Freedom  Accelero¬ 
meters  with  Mutually  Perpendicular  Input  Axes. 


sultant  specific  force  on  the 
plane  normal  to  the  axis  of  free¬ 
dom  of  the  pendulum.  The  out¬ 
put  of  such  a  unit  is  the  angle 
between  the  indicated  vertical 
and  the  pendulum  tracking  line 
about  this  axis.  Two  such  units, 
with  their  axes  of  freedom  non¬ 
parallel  (preferably  mutually 
perpendicular)  in  the  indicated 
horizontal  plane,  are  required 
for  the  complete  indication  of 
the  vertical. 

For  comparison,  a  vector 
diagram  associated  with  the 
operation  of  a  linear  accelero¬ 
meter  is  given  in  Fig.  6.  A 
single-degree-of-freedom  ac¬ 
celerometer  measures  the  com¬ 
ponent  of  the  resultant  specific 
force  along  its  own  input  axis. 
The  output  of  the  unit  is  a  signal 
proportional  to  the  magnitude 
of  this  component.  As  with 
single-degree-of-freedom  pen¬ 
dulums,  two  such  units,  with 
their  input  axes  preferably  mu¬ 
tually  perpendicular  in  the  in¬ 
dicated  horizontal  plane,  are 
required  for  the  full  indication 


of  the  vertical.  (Various  methods  of  controlling  the  orientation  of  the  plat¬ 


form  on  which  the  accelerometer  or  pendulum  are  mounted,  on  the  basis  of 
signals  associated  with  the  output  of  either  device,  are  subsequently  dis¬ 
cussed.) 


Note  that  the  information  obtainable  from  either  pendulums  or  accel 
erometers  is  essentially  the  same.  The  nature  of  the  information  is 


*  The  concept  of  an  orientation  as  the  mechanical  input  or  output  of  a  component  deserves  am¬ 
plification  here,  since  it  will  be  applied  repeatedly  in  this  report.  An  orientation  is  defined  as  an 
angular  displacement,  and  in  the  functional  diagrams  of  this  report  includes  all  time  derivatives 
of  the  displacement;  i.e.,  instantaneous  angular  velocity,  angular  acceleration,  etc. 


Figure  7.  Line  Schematic  Diagram  Showing  Essential  Features  of  a  Single-Refer¬ 
ence-Direction  Stabilization  Control  Unit  Using  a  Two-Degree-of-Freedom  Gyro¬ 
scopic  Element. 
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different:  a  pendulum  tracks  a  direction  and  an  accelerometer  measures  a 
magnitude.  Still,  no  basic  information  is  obtainable  from  the  use  of  one 
type  of  unit  that  cannot  also  be  obtained  from  the  use  of  the  other.  For  this 
reason,  it  is  not  possible  to  use  accelerometer  data  to  correct  pendulum 
data,  or  vice  versa.  The  choice  of  unit  is  dictated  by  engineering  design 
and  operating  requirements  rather  than  by  any  basic  difference  in  the  data 
obtained. 

It  has  been  shown(3)  that  vertical-indicating  performance  is  equivalent 
whether  pendulums  or  accelerometers  are  used  to  obtain  specific  force 
data.  With  an  accelerometer  the  Schuler  tuning  condition  depends  only  on 
the  magnitude  of  the  radius  of  the  Earth,  and  the  system  dynamics  (natural 
period,  damping,  etc.)  depend  on  the  magnitude  of  the  vertical  component 
of  the  resultant  specific  force.  With  a  pendulum,  the  Schuler  tuning  condition 
depends  on  the  magnitudes  of  both  the  radius  of  the  Earth  and  the  vertical 
component  of  the  resultant  specific  force,  but  the  system  dynamics  are  in¬ 
dependent  of  either  of  these  two  factors.  In  all  other  performance  character¬ 
istics  the  vertical-indicating  systems  are  similar. 

Gyroscopes;  Basic  Considerations 

A  detailed  account  of  the  performance  of  a  gyroscope  is  given  in  a  paper 
by  J.  J.  Jarosh(26)_  The  essential  geometrical  features  associated  with  gyro¬ 
scope  operation  are  illustrated  in  Figs.  7  and  8.  A  two-degree-of-freedom 
gyroscope  (Fig.  7),  to  which  no  torques  are  deliberately  applied  and  which 
has  negligibly  small  uncertainty  torques,  will  maintain  the  direction  of  the 
spin  axis  of  its  rotor  fixed  with  respect  to  inertial  space.  *  This  is  merely 
a  manifestation  of  Newton's  first  law  of  motion  applied  to  rotation.  In  this 
case,  the  angular  momentum  of  the  gyro  directly  furnishes  the  torque  against 
which  orientation  stabilization  can  be  made.  The  spin  axis  of  such  a  gyro 
unit  then  becomes  equivalent  to  a  stellar  line  of  sight,  a  property  that  is  of 
prime  importance  to  inertial  navigation.  When  a  torque  is  applied  to  a 
gyroscope,  the  gyro  rotor  unit  rotates  so  as  to  align  its  spin  axis  with  the 
applied  torque  at  a  rate  proportional  to  the  magnitude  of  the  torque.  This 
rotation  is  called  precession,  and  is  a  manifestation  of  Newton's  second 
law  of  motion  applied  to  rotation.  The  gyro  unit  under  such  conditions 
serves  as  an  integrator,  yielding  a  rate  of  change  of  angular  position  as 
output  with  a  torque  as  input.  Follow-up  signal  generators  are  needed  to 
furnish  information  on  the  orientation  of  the  controlled  platform  that  is 
using  the  two-degree-of-freedom  gyro  unit  as  its  angular  orientation 
reference. 

A  single -degree -of -freedom  gyroscope  must  operate  in  conjunction 
with  a  closed -loop  power -driven  platform  (Figure  8)  that  supports  the 
gyro  and  is  oriented  on  the  basis  of  data  received  from  the  gyro.  A  care¬ 
ful  distinction  must  be  drawn  between  this  gyro  application  and  the  one 


*  Actually,  it  is  the  angular  momentum  that  remains  fixed  in  snace.  In  a  practical  gyro  unit,  the 
mass  distribution  is  such  that  the  spin  axis  of  the  rotor  is  effectively  a  principal  axis  of  inertia, 
and  nutation  is  negligible;  so  that  the  angular  momentum  vector  and  the  spin  axis  of  rotor  are  effec¬ 
tively  parallel. 
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just  discussed,  because  here  the  gyro  units  themselves  do  not  exert  effective 
stabilizing  torques;  they  merely  act  as  signal  generators.  The  performance 
of  such  a  stabilization  system  depends  on  whether  the  gyroscope  is  of  the  inte¬ 
grating  type  or  the  rate  type.  An  integrating  gyro  to  which  no  torques  are 
deliberately  applied*,  and  which  has  negligibly  small  uncertainty  torques,  will 
maintain  by  means  of  its  attendant  closed-loop  system  an  arbitrary  direction 
fixed  with  respect  to  inertial  space,  in  the  plane  normal  to  the  input  axis** 
of  the  gyro  unit.  Two  such  units,  with  their  input  axes  preferably  mutually 
perpendicular,  can  cause  the  direction  of  the  normal  to  the  plane  of  their 
input  axes  to  remain  fixed  with  respect  to  inertial  space.  This  means  that 
two  single-degree-of -freedom  integrating  gyro  units  can  cause  a  direction  t,o 
remain  fixed  in  inertial  space  in  a  manner  equivalent  to  the  maintenance  of 
the  direction  by  a  single  two-degree-of -freedom  gyro  unit. 

A  single-degree-of -freedom  rate  gyro  operates  in  a  manner  similar  to  the 
analogous  integrating  gyro  except  that  the  maintenance  of  a  direction  in  space 
depends  on  the  nulling  of  the  angular  velocities  acting  on  the  platform  rather 
than  on  the  nulling  of  angular  displacements.  This  mode  of  operation  can  lead 
to  a  drift  of  the  direction  that  is  being  maintained  away  from  the  reference  dir¬ 
ection  fixed  in  space.  A  platform  controlled  by  single-degree-of -freedom  gyro 
units  (either  rate  or  integrating;  see  Fig.  8)  will  rotate  at  a  rate  proportional 
to  the  torque  applied  to  the  gyro  units,  thus  furnishing  integrator  or  velocity- 
drive  characteristics. 

The  method  of  using  a  gyro  unit  for  space  stabilization  data  has  a  marked 
effect  on  the  theory  of  operation  of  the  device.  In  one  class  of  systems  the 
gyros  are  left  essentially  undisturbed  and  maintain  a  reference  direction  in 
inertial  space.  In  the  other  class  of  systems,  torques  are  applied  to  the  gyros 
to  obtain  a  direction  controlled  by  data  from  the  system's  tracking  unit. 

Methods  for  Indicating  Astronomical  Position 

Basic  Properties  of  a  Feasible  Vertical  Indicator 

In  the  preceding  discussion  of  the  Schuler  pendulum,  it  was  brought  out  that 
it  is  not  physically  feasible  to  use  a  pendulum  directly  as  a  precise  indicator  of 
the  vertical  under  dynamic  conditions.  A  practical  system  for  indicating  the 
vertical  under  dynamic  conditions  could,  however,  be  made  from  a'  closed  servo 
loop  incorporating  a  form  of  Schuler  tuning  in  the  loop  dynamics.  The  essential 
geometrical  parts  of  such  a  loop  are  shown  in  Figure  9.  This  device  will  be 
called  a  vertical  indicator.  Its  function  is  to  maintain  a  platform  called  the 
controlled  member  in  the  horizontal  plane.  No  further  stabilization  will  be 
discussed  here;  specifically,  the  subject  of  the  compass-like  horizontal 
orientation  of  the  platform  will  be  deferred  for  later  examination.  There 
are  three  interrelated  subsystems  in  the  vertical  indicator,  two  of  which  are 


*  Torque  supplied  by  rotation  of  the  platform  supporting  the  gyro  unit  are  not  classed  here  as 
deliberately  applied  torques;  they  are  inherent  results  of  the  stabilizing  process  and  cause  null  po¬ 
sition  oneration  of  the  gyro  unit. 

**  The  input  axis  of  a  single-degree-of-freedom  gyro  unit  is  the  axis  about  which  an  angular  velo¬ 
city  of  the  base  supporting  the  gyro  unit  can  cause  a  precession  of  the  gyro  spin  axis  about  its  single 
axis  of  freedom. 
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displayed  in  Figure  9.  The  three  subsystems  are  as  follows: 

1.  a  system  for  obtaining  resultant  specific  force  data,  involving  either 
accelerometers  or  pendulums; 

2.  a  system  for  orienting  a  controlled  member  in  response  to  data  from 
(1)  to  provide  the  vertical  indication,  involving  gyro  units*  and  servo 
drives;  and 

3.  a  system  for  properly  modifying  and  filtering  data  from  (1)  above 
before  they  are  applied  to  (2)  in  order  to  control  the  dynamic  per¬ 
formance  of  the  whole  system. 

The  vertical  indicating  system  shown  schematically  in  Fig.  9  shows 
accelerometers  and  single-degree-of-freedom  gyro  units  for  subsystems  (1) 
and  (2).  Subsystem  (3)  would  consist  of  associated  electronic  and  electro¬ 
mechanical  equipment  not  shown  in  Fig.  9. 

The  Integration  of  Rates 

When  a  system  of  the  kind  shown  in  Fig.  9  is  carried  over  the  Earth’s  sur¬ 
face,  it  follows  a  path  that  involves  curvature  with  respect  to  space.  In  the  case 
of  the  ideally  spherical  Earth,  the  path  can  be  a  great  circle.  A  simple  picture 
of  the  role  rate  integration  plays  in  the  indication  of  astronomical  position  is 
obtained  by  supposing  the  vehicle  to  have  a  meridian  course  (Fig.  10  ).  The 
controlled  member  of  Fig.  9  is  started  at  the  equator,  the  zero  reference  for 
latitude^,  as  the  controlled  member  accelerates  northward  along  a  meridian  in 
the  plane  of  Fig.  10,  it  suffers  an  essentially  geocentric  angular  acceleration. 
Thus  the  controlled  member  itself  must  be  re -oriented  about  its  axes  with  re¬ 
spect  to  the  Earth,  For  simplicity,  suppose  that  its  x-axis  (Figure  9)  is  in  the 
plane  of  Fig.  10;  the  controlled  member  then  has  an  angular  acceleration 
about  its  y-axis,  which  is  perpendicular  to  the  plane  of  Fig.  10.  It  receives 
this  angular  acceleration  essentially  in  response  to  the  sensing  element  moun¬ 
ted  on  it  —  an  accelerometer  or  pendulum.  Either  device  responds  to  the  x-com 
ponent  of  the  inertia  reaction  specific  force  of  the  controlled  member  as  well 
as  being  simultaneously  responsive  to  the  x-axis  component  of  gravity. 

A  simplified  vertical  indicator  which  would  serve  this  particular  purpose 
is  shown  in  Figure  11.  The  accelerometer  necessarily  detects  linear  acceler¬ 
ations,  but  the  information  from  it  must  be  processed  to  obtain  an  essentially 
geocentric  angular  displacement.  This  processing  therefore  involves  a  double 
integration  with  respect  to  time.  The  fact  that  the  incremental  angular  dis¬ 
placement  is  essentially  geocentric  implies  that  a  knowledge  of  the  average 
Earth  radius  over  the  course  would  be  required  to  furnish  the  distance  of  travel 
of  the  vehicle  on  the  Earth,  i.  e. ,  one  minute  of  arc  corresponds  to  one 
nautical  mile. 


*  In  an  actual  system,  the  gyro  units  are  part  of  the  means  for  isolating  the  controlled  mem¬ 
ber  from  arbitrary  rotations  of  the  supporting  structure  and  vehicle,  i.  e.,  roll  and  pitch;  this  is 
taken  up  subsequently.  Meanwhile,  the  role  of  the  gyro  units  in  the  systems  about  to  be  discussed 
will  be  confined  to  their  function  as  part  of  the  controlled  member  drive  system. 
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The  integration  of  the  accelerometer  output  data  can  be  accomplished 
in  two  fundamentally  different  ways.  These  will  be  called  analytical 
integration  and  geometric  integration,  respectively.  Analytical  integration 
is  characterized  by  modifying  the  accelerometer  output  data  by  means  of 
two  integrators,  so  that  the  geocentric  angular  displacement  is  simply  the 
output  of  the  second  integrator.  The  first  integrator  receives  the  angular 
acceleration  of  the  indicated  vertical  as  an  input  and  delivers  the  angular 
velocity  of  the  indicated  vertical  to  the  second  integrator.  As  was  men¬ 
tioned  previously  in  connection  with  the  role  of  vertical  indication  in 
navigation  systems,  the  accelerometer  unit  must  be  mounted  on  a  stable 
platform,  with  its  input  axis  tangential  to  the  Earth,  or  it  will  indicate 
gravity  components  as  well  as  the  vehicle's  linear  acceleration  on  its 
course.  This  stabilization  of  the  accelerometer  mounting  may  then  be 
regarded  as  a  correction  to  the  accelerometer  input,  since  if  the  platform 
horizontal  stabilization  system  is  Schuler-tuned,  it  will  tend  to  remove 
gravity  effects  from  the  accelerometer  input. 


Geometric  integration,  on  the  other  hand,  is  characterized  by  comparing 
the  indicated  vertical  with  a  reference  vertical.  In  the  simple  case  under 
examination,  the  reference  vertical  could  be  the  direction  of  the  indicated 
vertical  at  the  equator.  If%his  reference  is  physically  carried  with  the  sys¬ 
tem  as  it  moves  north,  latitude  will  be  indicated  directly  as  the  angle  be¬ 
tween  the  indicated  vertical  and  the  reference  vertical.  Evidently,  the  in¬ 
dicated  vertical  must  be  extracted  as  a  datum  from  this  system  in  order  that 
the  system  indicate  latitude. 

Other  essential  differences  between  analytical  and  geometric  integration 
of  the  accelerometer  output  data  will  be  examined  further  in  the  discussion 
following. 
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The  Basic  Vertical  Indicator 

The  method  of  integration  for  position  information  is  immaterial  as  far 
as  the  orientation  of  the  controlled  member  mounting  the  accelerometer  is 
concerned:  the  controlled  member  must  indicate  the  vertical*.  Vertical  in¬ 
dication  is  fundamental  to  the  indication  of  geocentric  angular  displacement; 
the  nlatform  mounting  the  accelerometer  must  be  tangential  to  the  Earth. 

Before  discussing  various  methods  for  determining  the  geocentric  angular 
displacement,  therefore,  it  will  be  of  value  to  examine  in  some  detail  a  sim¬ 
plified  (i.e.,  ideal,  undamped)  single-degree-of-freedom  vertical  indicator 
of  the  kind  sketched  in  Fig.  11.  The  extension  of  the  following  arguments  to 
the  two-degree-of -freedom  vertical  indicator  (Fig.  9)  then  follows  immed¬ 
iately.  A  functional  diagram  for  the  essential  operating  components  of  the 
single-degree-of-freedom  system  of  Fig.  11  is  shown  in  Fig.  12.  This  dia¬ 
gram  is  essentially  duplicated  as  Fig.  3-2  of  Derivation  3  appended.  For  the 
purpose  of  mathematical  analysis  there,  a  symbolic  representation  of  the 
performance  characteristics  of  the  components  is  used  to  label  them.  In 
Derivation  3  the  performance  equations  for  this  system  are  derived  from 
basic  physical  considerations.  The  reader  is  referred  there  for  a  mathe¬ 
matical  substantiation  of  the  conclusions  drawn  in  connection  with  the  present 
discussion.  It  should  be  noted  that  only  ideal,  undamped  systems  are  dis¬ 
cussed  in  this  report;  non-ideal,  damped  systems  are  considered  in  the  next 
report. 

Schuler  Tuning  in  a  Closed-Loop  System 

The  performance  equation  for  the  system  diagrammed  in  Figure  12  is  ex¬ 
pressed  in  Derivation  3  in  terms  of  the  correction  to  the  indicated  vertical, 
which  is  the  angle  between  the  indicated  vertical  and  the  true  vertical,  i.e. , 
this  correction  is  the  negative  of  the  inaccuracy  in  indicating  the  vertical.  The 
performance  equation  is  a  second-order  differential  equation  which  will  con¬ 
tain  a  term  in  the  angular  acceleration  of  the  true  vertical  unless  the  coefficient 
of  this  term  is  equal  to  zero.  This  term  is  part  of  the  description  of  a  system 
in  which  the  indicated  vertical  consistently  lags  or  leads  the  true  vertical,  and 
equating  the  coefficient  to  zero  tends  to  remove  the  lead  or  lag,  thus  accomplish¬ 
ing  an  effective  Schuler  tuning  of  the  system.  Equating  the  coefficient  of  the 
true  vertical's  angular  acceleration  to  zero  amounts  to  adjusting  the  sensitivities 
of  three  components  of  the  system  so  that  their  product  is  the  reciprocal  of  the 
average  Earth  radius.  This  adjustment  is  applied  to  each  of  the  systems  to  be 
discussed  to  obtain  Schuler  tuning. 

In  Derivation  3,  an  equation  is  derived  which  expresses  the  cori'ection  to  the 
indicated  vertical  as  a  function  of  the  time.  As  would  be  expected,  considering 
the  assumption  that  the  system  is  undamped  and  Schuler-tuned,  this  performance 
equation  is  essentially  that  for  a  simple  harmonic  oscillator,  with  a  period  equal 
to  the  Schuler  period,  and  with  displacement  and  velocity  amplitudes  controlled 
by  the  initial  conditions  of  the  motion  —  i.e. ,  the  coefficient  of  the  cosine  term 

*  Except  for  certain  short  duration  situations  not  applicable  to  general  submarine  navigation. 
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is  the  correction  to  the  indicated  vertical  at  tire  start  of  the  problem,  and  the 
coefficient  of  the  sine  term  involves  the  time  rate  of  change  of  this  correction 
at  the  start  of  the  problem. 

It  should  be  noted  here  that  a  superficially  different  point  of  view  is  pos¬ 
sible  concerning  the  operation  of  a  vertical  indicator:  it  may  be  said  to  "track" 
the  indicated  vertical.  Resultant  acceleration  is  the  basic  input  to  an  acceler¬ 
ometer  or  a  pendulum.  Now  when  a  pendulum  support  is  accelerated,  -the 
pendulum  bob  tends  to  lag  behind  the  support.  In  the  most  general  case,  this 
"lag"  of  the  pendulum  bob  results  in  indication  of  the  apparent  vertical  instead 
of  the  true  vertical,  where  the  apparent  vertical  is  the  direction  of  the  resultant 
specific  force  and  lags  the  true  vertical  by  an  angle  that  is  a  function  of  the 
acceleration  acting.  Consider  that  the  vertical  indicating  system  is  a  tracking 
device  that  receives  resultant  specific  force  as  its  input  and  gives  the  indicated 
vertical  as  its  output.  The  input  receiving  member  is  a  damped  pendulum, 
which  is  a  very  satisfactory  "tracker"  of  the  apparent  vertical.  The  controlled 
member  is  the  tracking  member;  the  direction  of  a  normal  to  the  plane  estab¬ 
lished  by  the  pendulum  unit  input  axis  is  defined  as  the  indicated  vertical.  If 
the  indicated  vertical  is  to  be  parallel  to  the  true  vertical,  it  is  essential  that 
the  indicated  vertical  be  displaced  from  the  apparent  vertical  by  some  angle. 

The  same  considerations  apply  if  an  accelerometer  replaces  the  pendulum 
mentioned  above.  From  this  point  of  view,  the  vertical  indicator  is  a  closed- 
loop  system  with  a  "lag"  between  output  and  input,  which  is  a  fundamental  and 
necessary  characteristic  of  its  operation. 

Fundamentals  of  Position  Indication:  Summary 

Astronomical  position  measurement  involves  the  determination  of  the 
essentially  geocentric  angles  between  indicated  and  reference  verticals.  It 
is  proposed  to  do  this  with  a  stabilized  platform  as  an  essential  feature.  The 
stabilization  loop  receives  specific  force  data  from  accelerometers  or  pen¬ 
dulums,  and  delivers  these  data  to  gyro  units  and  servo  drives,  which  in  turn 
drive  the  platform  toward  a  horizontal  position.  The  specific  force  data  are 
modified  and  filtered  to  control  the  system  dynamics.  Part  of  this  modifi¬ 
cation  consists  of  Schuler -tuning  the  stabilization  loop.  Two  basic  methods 
of  utilizing  the  indicated  vertical  are  then  available  to  obtain  the  vehicle's 
essentially  geocentric  angle  of  travel.  One,  called  here  analytical  integration, 
involves  the  integration  of  accelerometer  outputs  or  of  gyro  inputs,  with  the 
accelerometers  and  gyro  units  mounted  on  the  stable  platform.  The  other 
method,  called  here  geometric  integration,  involves  the  direct  comparison 
of  the  indicated  vertical  with  a  reference  vertical.  Both  methods  will  now 
be  outlined  in  terms  of  six  representative  systems. 


Typical  Methods  for  Inertial  Indication  of  Position 

Open-Chain  Integration  of  the  Angular  Velocity  of  the  Indicated  Vertical 

Figure  13  is  a  functional  diagram  of  the  essential  operating  components 
of  a  latitude  indicator  which  uses  an  open-chain  single  integration  of  the  an¬ 
gular  velocity  of  the  indicated  vertical  to  arrive  at  the  indicated  latitude. 
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Figure  13  corresponds  to  Fig.  4-1  in  Derivation  4.  Note  that  this  system  is 
primarily  the  basic  vertical  indicator  of  Fig.  12,  with  the  angular  velocity  in¬ 
tegration  tapped  off  the  vertical  indicator  loop  at  the  angular  velocity  signal 
that  calls  for  precessing  the  gyro  unit,  so  that  the  position  integrator  output  is 
the  indicated  change  in  latitude.  The  vertical  indicator  is  undisturbed  by  the 
presence  of  this  integrator,  as  far  as  the  loop  dynamics  are  concerned.  Hence, 
(as  is  shown  in  Derivation  4)  the  same  Schuler  tuning  condition  applies  to  this 
system  as  applied  to  the  basic  vertical  indicator  of  Fig.  12. 

A  two-degree-of -freedom  system  constructed  along  these  lines  would 
require  an  x-axis  accelerometer  in  addition  to  the  y-axis  accelerometer 
shown  in  Fig.  13,  as  well  as  x-axis  components  corresponding  to  the  orien¬ 
tational  control  signal  generating  system,  gyro  unit,  and  controlled  member 
drive.  In  addition,  the  x-axis  system  requires  that  its  angular  velocity  sig¬ 
nal  be  modified  by  secant  of  latitude  in  order  to  give  longitude  rate  data  to 
its  position  integrator.  The  position  integrator  in  the  x-axis  system  would 
indicate  longitude  difference  with  respect  to  some  arbitrary  reference.  The 
maintenance  of  the  y  and  x  axes  along  a  parallel  of  latitude  and  a  meridian, 
respectively,  must  be  handled  by  an  auxiliary  orientation  system,  e.g.,  a 
gyrocompass. 

The  gyro  units  ideally  maintain  their  spin  axes  fixed  with  respect  to  in¬ 
ertial  space.  To  maintain  the  controlled  member  on  which  they  are  mounted 
tangential  to  the  Earth  in  its  daily  rotation,  therefore,  the  x-axis  gyro  must 
be  precessed  by  a  correction  torque  which  keeps  the  platform  up  with  Earth 
rate.  The  calculations  of  Derivation  4  consider  this  torque  to  be  generated 
by  a  special  precession  current  applied  to  the  torque  generator  of  the  gyro 
unit. 


In  Derivation  4,  it  is  shown  that  if  the  controlled  member  has  the  equation 
of  motion  of  a  Schuler-tuned  pendulum, the  correction  to  the  indicated  position 
(that  is,  the  instantaneous  difference  between  the  true  and  indicated  positions) 
has  a  similar  equation  of  motion,  except  for  the  appearance  of  certain  initial 
conditions  as  additive  constants  of  integration  in  the  indicated  position.  Spec¬ 
ifically,  the  correction  to  the  indicated  position  at  any  time  is  the  correction 
to  the  indicated  position  at  the  start  of  the  problem,  plus  the  present  correc¬ 
tion  to  the  indicated  vertical,  less  the  initial  correction  to  the  indicated  ver¬ 
tical.  The  indicated  position  is  thus  Schuler -periodic  with  respect  to  boundary¬ 
matching  terms  associated  with  the  vertical  indication. 

Open-Chain  Double  Integration  of  the  Angular  Acceleration  of  the  Indicated  Vertical 

Figure  14  is  a  functional  diagram  of  the  essential  operating  components 
of  a  latitude  indicator  which  uses  an  open-chain  double  integration  of  the  an¬ 
gular  acceleration  of  the  indicated  vertical  to  arrive  at  the  indicated  latitude. 
Figure  14  corresponds  to  Fig.  5-1  of  Derivation  5.  This  system  represents 
only  a  slight  modification  of  that  of  Fig.  13,  and  the  prior  discussion  con- 
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the  Angular  Acceleration  of  the  indicated  Vertical. 
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cerning  Schuler  tuning,  the  extension  two  degrees  of  freedom  and  longitude  as 
well  as  latitude  indication,  and  the  gyrocompass  requirement,  is  applicable. 

It  should  be  noted  that  while  the  accelerometer-unit  itself  indicates  linear 
acceleration,  the  desired  system  output  is  an  essentially  geocentric  angle  of 
travel.  The  accelerometer  output  is,  in  fact,  proportional  to  an  angular  ac¬ 
celeration,  namely,  the  linear  acceleration  of  the  vehicle  divided  by  the  average 
Earth  radius  over  the  course.  This  Earth-radius  value  is  important  in  the 
Schuler  tuning  of  the  system,  as  described  in  Derivation  3,  and  determines 
the  dynamics  of  the  vertical  indicator.  However,  it  is  emphasized  that  the 
Earth  radius  has  no  other  role  in  the  indication  of  position.  In  effect,  by 
stabilizing  the  platform,  the  Schuler  tuning  condition  furnishes  the  average 
Earth  radius  at  which  the  accelerometer  unit  operates,  so  that  an  essentially 
geocentric  angular  acceleration  is  immediately  derivable  from  its  output. 

The  performance  equation  for  this  system  is  given  in  Derivation  5. 

The  correction  to  the  indicated  position  involves: 

(a)  the  initial  correction  to  the  indicated  position 

(b)  the  initial  time  derivative  of  the  correction  to  the  indicated  position; 
this  term  will  cause  the  correction  to  increase  directly  with  the  time 
of  operation,  since  it  represents  a  false  ground  speed  component 

(c)  the  initial  correction  to  the  indicated  vertical;  the  correction  is 
Schuler-periodic  with  respect  to  this  term 

(d)  the  initial  time  derivative  of  the  correction  to  the  indicated  vertical; 
this  term  appears  as  part  of  two  functions:  the  first  will  produce 
Schuler-periodic  oscillations,  and  the  second  will  cause  the  correc¬ 
tion  to  increase  directly  with  time  of  operation. 

Direct  Double  Integration  of  Acceleration 

Figure  15  is  a  functional  diagram  of  the  essential  operating  components  of 
a  latitude  indicator  which  uses  a  doubly -integrating  accelerometer  as  an  inher¬ 
ent  component  in  the  vertical-indicating  loop.  The  output  of  such  an  accelero¬ 
meter  is  immediately  proportional  to  the  change  in  indicated  latitude.  However, 
this  output  is  not  a  suitable  orientational  control  signal  for  the  gyro  unit,  which 
must  be  driven  by  a  signal  corresponding  to  the  angular  velocity  of  the  indicated 
vertical.  The  vertical  indicating  loop  for  this  system,'  as  distinguished  from 
the  three  systems  previously  taken  up,  therefore  contains  a  differentiator  be¬ 
tween  the  accelerometer  and  the  gyro  unit.  This  slightly  modifies  the  expres¬ 
sion  for  the  Schuler  tuning  condition  for  this  loop,  in  that  the  differentiator 
sensitivity  appears  instead  of  the  integrator  sensitivity  used  in  the  loops  of 
Figs.  12,  13,  and  14.  Otherwise,  the  Schuler  tuning  condition  is  unchanged. 

The  extension  of  this  description  of  a  single -axis  system  to  that  involving  two 
degrees  of  freedom  goes  as  in  the  case  of  the  system  of  Fig.  14.  A  gyrocom¬ 
pass  or  similar  azimuth  device  is  also  a  requirement  here.  Derivation  6 
shows  that  the  equation  of  motion  of  this  position  indicator  is  the  same  as  for 
the  system  utilizing  Open-Chain  Double  Integration.  The  correction  to  the  in¬ 
dicated  position  will  again  involve  the  same  four  quantities  resulting  in  terms 
with  Schuler-periodic  oscillations  plus  terms' that  build  up  with  time. 
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Electromagnetic  Composition  of  Tangential 

and  Vertical  Components  of  Celestial  Longitude  Rate 

If  the  system  of  Fig.  9  is  stabilized  in  azimuth,  by  means  to  be  taken  up 
subsequently,  two  components  of  celestial  longitude  rate  (longitude  rate  taken 
with  respect  to  inertial  space)  may  be  used  to  obtain  indicated  latitude  and 
indicated  longitude  difference.  Figure  16  is  a  partial  functional  diagram  of 
such  a  system.  The  signals  feeding  the  x-gyro  unit  (Fig.  12)  and  the  z-gyro 
unit  (Fig.  9)  are  respectively  proportional  to  these  celestial  rate  components; 
the  x -signal  is  proportional  to  a  quantity  given  by  the  sum  of  Earth-rate  and 
longitude  rate  multiplied  by  the  cosine  of  indicated  latitude,  and  the  z-signal 
is  proportional  to  the  same  quantity  multiplied  by  the  sine  of  the  indicated 
latitude.  These  statements  are  justified  in  Derivation  8,  and  the  extraction 
of  position  information  from  these  data  is  also  covered  there  in  detail.  In 
brief,  the  method  is  as  follows:  the  aforementioned  cosine  and  sine  com¬ 
ponent  signals  are  fed  to  a  trigonometric  resolver  (see  Fig.  8-1  in  Deriva¬ 
tion  8).  Under  these  input  conditions  the  resolver  acts  as  a  device  which 
composes  the  sine  and  cosine  component  inputs  into  the  magnitude  of  the 
vector  (Earth-rate  plus  longitude  rate),  and  the  resolver  rotor  is  automat¬ 
ically  rotated  by  an  associated  servo  to  the  angle  corresponding  to  the  sine 
and  cosine  terms;  indicated  latitude,  in  tins  case.  The  magnitude  of  the 
vector  must  be  integrated  to  obtain  longitude  information.  Furthermore, 
the  Earth-rate  term  must  be  removed,  either  before  or  after  the  integra¬ 
tion,  to  make  this  information  represent  longitude  difference.  The  inte¬ 
grated  longitude  rate,  i.e.,  longitude  difference,  is  then  referred  to  a  ref¬ 
erence  longitude,  so  that  the  second  output  of  the  system  is  indicated  lon¬ 
gitude.  The  mechanization  of  an  ideal  system  performing  these  functions 
is  detailed  in  Fig.  8-2  of  Derivation  8. 

Geometric  Integration  Using  Angular  Velocity 

of  the  Indicated  Vertical  with  Respect  to  Inertial  Space  t 

Another  approach  to  the  position-indicating  problem,  using  geometric 
integration,  is  shown  in  the  form  of  a  functional  schematic  diagram  in  Fig.  17, 
corresponding  to  Fig.  9-1  of  Derivation  9.  The  method  is  as  follows:  A  stable 
platform  is  set  up,  indicating  the  vertical  —  essentially  the  two -axis  version 
of  Fig.  9.  This  platform  is  then  oriented  by  a  gyrocompass,  by  a  method  to 
be  described,  establishing  a  geographic  north  line.  On  this  stable  platform,  a 
latitude  gimbal  is  erected.  This  gimbal  has  an  axis  of  rotation  in  the  horizon¬ 
tal  east-west  direction  (Fig.  18)  and  is  automatically  elevated  about  this  axis 
on  the  basis  of  angular  velocity  information  in  the  indicated  meridian  plane. 
This  orientation  is  accomplished  by  the  latitude  indication  loop  (Fig.  17).  In 
equilibrium,  the  latitude  gyro  unit,  which  is  the  angular  velocity  tracker,  has 
its  input  axis  perpendicular  to  the  indicated  Earth's  polar  axis.  When  the  la¬ 
titude  gimbal,  on  which  the  gyro  unit  is  mounted,  assumes  an  orientation  other 
than  this  equilibrium  orientation,  the  gyro  unit  (a  rate  gyro)  receives  an  angu¬ 
lar  velocity  component  input.  The  gyro  unit  then  generates  a  signal  which, 
after  integration,  is  used  to  rotate  the  latitude  gimbal  to  its  equilibrium  orien¬ 
tation,  thereby  causing  null  operation  of  the  rate  gyro.  Latitude  is  then  directly 
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indicated  (Fig.  18)  as  the  angle  between  the  latitude  gimbal  and  the  indicated 
vertical.  The  performance  equation  for  this  system  is  developed  in  Deriva¬ 
tion  9. 

Within  the  latitude  gimbal,  the  longitude  gimbal  (Fig.  17)  rotates  about 
the  indicated  Earth's  polar  axis  to  register  indicated  celestial  longitude  differ¬ 
ence,  or  celestial  longitude  with  respect  to  some  arbitrary  sidereal  reference. 
The  longitude  gimbal  supports  the  longitude  gyro  unit  (an  integrating  gyro). 
The  gyro  unit,  which  has  its  input  axis  parallel  to  the  indicated  Earth's  polar 
axis,  receives  indicated  celestial  longitude  rate  as  an  input.  The  gyro  unit 
then  generates  a  signal  which  is  used  to  rotate  the  longitude  gimbal  relative 
to  the  latitude  gimbal  at  the  negative  of  the  indicated  celestial  longitude  rate. 
In  other  words,  the  longitude  gimbal  tends  to  remain  fixed  in  inertial  space. 

In  the  position-indicating  system  just  discussed,  an  electromagnetic  vector 
comnositor  yielded  celestial  longitude  rate  as  an  output,  which  was  then  pro¬ 
cessed  to  obtain  indicated  longitude  difference,  by  the  subtraction  of  Earth- 
rate  and  subsequent  integration.  In  the  present  case,  Earth-rate  compen¬ 
sation  is  also  required  to  obtain  longitude  difference*,  but  the  required  inte¬ 
gration  is  in'  vent  in  the  operation  of  the  longitude  drive  system.  The  per¬ 
formance  eqir  ion  for  tuis  system  is  also  developed  in  Derivation  9.  It  is 
shown  there  mat  the  correction  to  the  indicated  latitude  consists  of  sinusoidal 
terms  involving  a  period  dependent  on  the  angular  momentum  of  the  latitude 
gyro,  certain  system  sensitivities,  celestial  longitude  rate,  and  a  term  invol¬ 
ving  the  acceleration  in  latitude. 

It  is  to  be  noted  that  this  entire  system  is  self -settling;  the  only  initial 
setting  required  is  the  longitude  reference,  which  can  be  conveniently  made 
at  the  longitude  indicating  system. 

Geometric  Integration  thing  a  Pre-Aligned  Inertial  Reference 

Figure  19  is  a  functional  diagram  of  the  essential  operating  components 
of  a  latitude  indicator  which  compares  the  indicated  vertical  with  a  pre- 
aiigm..'  refeience  vertical  parallel  to  the  Earth's  polar  axis.  The  reference 
vertical  is  maintained  by  a  pre -aligned  rate -integrating  gyro  unit  and  its 
associated  servo  drive.  The  gyro  unit  is  left  undisturbed  during  the  marine 
operations,  thereby  maintaining  an  inertial  space  reference  orientation. 
Figure  19  corresponds  to  Fig.  10-1  in  Derivation  10.  The  superposition  of 
the  outputs  of  two  indicators  —  one  vertical  indicator  of  the  tracking  type 
shown  ir.  Fig.  12  and  one  fixed-axis  reference  indicator  maintained  by  an 
inertial  gyro  —  yields  indicated  co -latitude.  The  fixed -avs  indicator  fur¬ 
nishes  essentially  the  indicated  polar  axis. 

*  Two  practical  general  methods  of  compensation  are: 

1.  A  sidereal  time  drive  (clock),  used  to  rotate  the  longitude  gimbal  about  the  indicated  polar 
axis,  which  otherwise  indicates  celestial  longitude  difference,  to  remove  integrated  Earth-rate  and 
leave  longitude  difference, 

2.  An  angular  data  transmitter  on  the  longitude  gimbal,  .jsed  to  remove  celestial  longitude  data 
to  a  remote  point,  where  the  orientation  of  the  data  receiver  and  the  orientation  of  a  sidereal  clock 
shaft  feed  a  differential  whose  output  is  indicated  longitude  difference. 
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The  indicated  vertical  of  such  a  single-axis  system  is  given  by  the  orien¬ 
tation  of  the  outer  gimbal  in  a  two-gimbal  system.  The  accelerometer  is 
mounted  on  this  gimbal.  The  gyro  is  mounted  on  the  inner  gimbal,  where  it 
maintains  the  reference  vertical.  Thus  a  reference  vertical  is  physically 
"carried  along"  with  the  system  as  it  moves  over  the  Earth.  The  angle  be¬ 
tween  the  reference  vertical  and  the  indicated  vertical  is  the  indicated  co¬ 
latitude.  The  reference  gyro  must  be  pre-aligned  from  external  information 
with  its  input  axis  in  the  east  -  west  direction  and  its  output  axis  along  the 
polar  axis  at  the  start  of  operations.  The  presence  of  an  external  pre¬ 
aligned  reference  obviates  the  need  for  a  gyrocompass  as  an  integral  operat¬ 
ing  component  in  this  system. 

The  extension  of  this  system  to  the  practical  case,  using  three -axis 
stabilization,  is  not  as  simple  as  in  the  systems  previously  discussed.  Three 
undisturbed  (inertial)  gyros  are  then  required  on  the  innermost  gimbal,  their 
input  axes  corresponding  to  mutually  perpendicular  x,  y,  and  z  axes  respectively, 
with  the  z-axis  oriented  along  the  Earth's  polar  axis  and  the  x  and  y  axes  arbi¬ 
trarily  oriented,  but  fixed  in  inertial  space  in  a  plane  parallel  to  the  equatorial 
plane  (Fig.  20).  The  gimbal  is  then  driven  about  the  z-axis  with  a  sidereal 
(clock)  drive  that  removes  the  effect  of  the  Earth 's  sidereal  rotation  from  the 
inner  gimbal.  When  the  system  changes  longitude,  indicated  longitude  rate  is 
an  input  to  the  z-gyro  unit.  The  z-gyro  unit  then  generates  a  signal  which  pro¬ 
duces  an  additional  relative  rotation  between  the  inner  gimbal  and  its  support, 
at  the  negative  of  the  indicated  longitude  rate.  This  causes  the  gyro  input  axes 
in  the  indicated  equatorial  plane  to  remain  fixed  in  inertial  space;  i.e.,  they 
maintain  equatorial  "star  lines". 

A  total  of  five  gimbals  is  required  for  complete  base  motion  isolation. 
Starting  with  the  innermost  gimbal,  these  are: 

1.  The  inertial  space  gimbal,  mounting  the  three  inertial  gyros. 

2.  The  meridian  gimbal,  which  supports  the  inertial  space  gimbal  through 
a  differential,  so  that  the  meridian  gimbal  can  effectively  be  driven 
around  the  inertial  space  gimbal. 

3.  The  azimuth  gimbal,  which  is  the  controlled  member,  supporting 
two  accelerometer  units.  Each  accelerometer  unit  is  a  component 
in  a  loop  similar  to  Fig.  19. 

4.  The  pitch  gimbal,  which  supports  the  azimuth  gimbal  and  isolates 
the  system  from  the  pitch  motion  of  the  ship. 

5.  The  roll  gimbal,  which  supports  the  pitch  gimbal  and  isolates  the 

system  from  ship's  roll.  . 

The  angle  between  the  planes  of  the  meridian  and  azimuth  gimbals  gives 
the  indicated  co-latitude.  The  angle  between  the  inertial-space  and  meridian 
gimbals  gives  the  indicated  longitude  difference.  The  performance  equation 
for  a  single-axis  version  of  this  system  is  given  in  Derivation  10,  where  it 
is  shown  that  the  position -indicating  ability  of  an  ideal  undamped  system  of 
this  type  is  affected  by  the  accuracy  of  the  initial  gyro-unit  alignment. 
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Indication  of  Roll  and  Pitch 

The  functional  diagram  shown  in  Fig.  8-1,  Derivation  8,  shows  a  complete 
two-degree-of -freedom  vertical  indicator.  This  involves  some  components  not 
yet  discussed  but  nevertheless  fundamental  to  all  of  the  systems  described  — 
those  using  analytical  as  well  as  geometric  integration.  These  components  are 
associated  with  the  mounting  of  the  controlled  member  in  each  case  and  with  the 
way  the  ship's  motion  influences  the  system.  The  latter  concerns  the  indication 
of  roll,  pitch,  and  yaw;  only  roll  and  pitch  will  be  taken  up  here,  and  yaw  will  be 
discussed  with  the  gyrocompass  later. 

In  Fig.  8-1,  the  ship  orientation  enters  the  gimbal  support.  Inside  this 
support  are,  in  order  of  decreasing  size,  the  roll  gimbal,  the  pitch  gimbal,  and 
the  azimuth  gimbal.  The  azimuth  gimbal  corresponds  to  the  stable  platform 
previously  discussed.  The  fact  that  the  accelerometers  and  tracking  gyro  units 
are  mounted  on  the  azimuth  gimbal,  although  they  generate  signals  which  drive 
the  roll  and  pitch  gimbals,  necessitates  the  resolvers  and  the  so-called  secant 
computer  shown  in  Fig.  8-1.  It  should  be  emphasized  that  this  geometric  roll 
and  pitch  indication  is  applicable  to  all  of  the  systems  already  discussed. 

The  Gyrocompass 

The  note  in  the  lower  left  corner  of  Fig.  9  states  that  "for  three -degree - 
of -freedom  stabilization  and  orientational  control,  the  z-axis  gyro  unit  must 
be  supplied  a  proper  orientational  control  signal  from  a  source  not  shown  in 
this  diagram".  This  source  will  necessarily  be  a  compass;  for  the  present 
purpose,  a  gyrocompass  is  indicated.  The  problem  in  orientation  here  is 
considerably  simplified  once  the  vertical-indication  problem  has  been  solved 
independently.  Standard  gyrocompass  methods  will,  in  general,  be  less  use¬ 
ful  than  the  one  to  be  described,  which  takes  full  advantage  of  the  horizontal 
stability  of  the  controlled  member  which  is  to  be  stabilized  in  azimuth. 

The  Earth  makes  one  revolution  about  its  axis  with  respect  to  the  "fixed 
stars"  every  sidereal  day  (23  hours,  56  minutes,  4.1  seconds  of  standard 
solar  time).  Hence,  it  possesses  an  angular  velocity  with  respect  to  inertial 
space.  At  any  point  on  the  surface  of  the  Earth,  this  angular  velocity  has 
components  which  are  directed  along  the  local  vertical  and  in  the  horizontal 
plane  directed  north.  Figure  21  shows  these  components  as  they  exist  at 
any  point  on  the  Earth's  surface.  Consider  the  special  case  where  the  sys¬ 
tem  of  Fig.  9  is  at  rest  with  respect  to  the  Earth,  and  is  settled,  as  far  as 
its  vertical  -  indicating  function  is  concerned.  Supoose,  however,  that  the 
input  axis  of  the  y-gyro  unit  is  not  pointed  east,  so  that  the  x-axis  reference 
line  on  the  controlled  member  does  not,  in  consequence,  indicate  true  north. 

The  angle  between  indicated  and  true  north  can  then  be  defined  as  the  correc¬ 
tion  to  indicated  north.  In  operation,  this  angle  will  always  be  small,  since 
the  function  of  the  azimuth  tracking  system  is  to  cause  the  controlled  member 
to  align  its  x-axis  with  true  north.  When  the  system  is  at  rest  on  the  Earth 
and  has  an  azimuth  error,  there  is  an  angular  velocity  component  input  to  the 
y-gyro  unit.  This  input  is  specifically  the  horizontal  component  of  Earth- 
rate,  projected  on  the  y-axis,  i.e,,  multiplied  by  the  sine  of  the  correction 
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Figure  20.  Diagram  Showing  Essential  Components  of  a  Three-Axis  Position  Indica¬ 
tor  Using  an  Inertial  Reference. 
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Figure  21.  The  Earth's  Angular  Velocity  Components  Which  Affect  a  Gyrocompass. 


to  indicated  north  (Fig.  22).  This  input  affects  the  y-axis  base  motion  isola¬ 
tion  system  (Fig.  12)  in  such  a  manner  as  to  cause  the  x-axis  to  dip  when 
this  axis  is  pointing -west  of  north,  i.e.,  the  y-axis  system  is  attempting  to 
stabilize  the  controlled  member  in  inertial  space.  This  dip  action  causes  the 
y-axis  pendulum  or  accelerometer  unit  to  send  an  orientational  control  signal 
to  the  y-gyro  unit,  thereby  precessing  the  controlled  member  relative  to  iner¬ 
tial  space  in  order  to  keep  the  controlled  member  horizontal.  Since  this  orien- 
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Figure  22.  Horizontal  Earth-Rate  Component  Projection  on  Y-Axis  Due  to  Azimuth  Misalignment. 


tational  control  signal  is  proportional  to  the  angular  velocity  projection  on  the  y 
axis,  it  is  a  measure  of  the  azimuth  error. 

Figure  23  shows  this  signal  as  an  input  to  the  z-axis  orientational  control 
signal  generating  system.  The  z-axis  controlled  member  drive  system,  on  the 
basis  of  this  Earth-rate  signal  from  the  y-gyro  unit,  then  orients  the  controlled 
member  in  azimuth  so  that  the  correction  to  indicated  north  tends  toward  zero; 
i.e.,  so  that  the  y-gyro  unit  Earth-rate  component  input  is  reduced  to  zero. 
However,  as  noted,  the  vertical  indicator  is  simultaneously  influenced  by 
this  input.  The  y-axis  portion  of  the  vertical  indicator  and  the  z-axis  sys¬ 
tem  are  coupled  in  such  a  way  that  their  dynamic  stability  characteristics  are 
mutually  interdependent.  It  is  shown  in  Derivation  7  that  their  performance 
equations  are  similar,  and  that  an  azimuth  disturbance  will  create  required 
corrections  to  the  indicated  vertical  and  to  indicated  north  which  are  approx¬ 
imately  Schuler  periodic. 


It  should  be  noted  that  a  separate  rate  gyro  unit,  with  its  input  axis  along 
the  y-axis,  rather  than  the  y-gyro  unit  of  Fig.  9,  would  accomplish  the  same 
purpose.  However,  the  additional  gyro  unit  would  be  unnecessary,  since  the 
y-gyro  unit  is  already  required  for  vertical  indication,  and  the  use  of  a  sep¬ 
arate  rate  gyro  unit  would  still  impose  simultaneous  stability  conditions  on  the 
y-  and  z-loops. 

Suppose  that  the  system  of  Fig.  9  moves  with  respect  to  the  Earth,  and 
that  the  z-axis  gyro  unit  of  Fig.  9  is  a  component  of  the  azimuth  stabilizing 
loop  of  Fig.  23.  The  latitude  rate  of  the  system  now  appears  in  the  y-axis 
orientational  control  signal  as  a  term  that  is  undesired  in  the  z-axis  system. 
Latitude  rate  can  be  effectively  removed  by  the  introduction  of  the  ship's 
northward  velocity  component  (from  an  external  source)  into  the  z-axis  orien¬ 
tational  control  signal  generating  system  (Fig.  23).  (This  correction  is  similar 
to  those  required  in  standard  gyrocompasses.) 

Change  of  longitude  involves  an  associated  rotation  of  the  meridian  about 
the  vertical  axis,  due  to  the  convergence  of  the  meridian  grid  poleward  (this 
effect  is  zero  at  the  equator).  Longitude  rate  effects*can  be  compensated  by 
injecting  the  x-axis  orientational  control  signal,  modified  by  the  tangent  of 
latitude,  into  the  z-axis  orientational  control  signal  generating  system.  If 
the  x-signal  has  in  addition  been  compensated  independently  for  Earth-rate, 
the  z-gyro  unit,  operating  with  respect  to  inertial  space,  will  also  be  so  com¬ 
pensated.  The  role  of  the  z-axis  orientational  control  signal  generating  sys¬ 
tem  in  further  processing  of  the  y-axis  orientational  control  signal  is  detailed 
in  Derivation  7. 

The  combined  system  described  thus  performs  the  dual  function  of  track¬ 
ing  the  true  vertical  and  true  north. 


SUMMARY 


The  self-contained  inertial  navigation  systems  discussed  in  this  report 
have  been  considered  only  in  ideal  form,  on  the  basis  of  the  theoretical  back¬ 
ground.  Imperfections  due  to  actual  equipment  are  deferred  for  the  next  re¬ 
port.  Accordingly,  no  attempt  is  made  here  to  find  a  most  feasible  system; 
this  also  is  left  to  the  next  report. 


APPENDIX 


The  mathematical  derivations  on  the  following  pages  are  referred  to  by 
number  in  the  foregoing  text.  The  self-defining  notation  used  in  the  equations 
is  based  on  a  formulation  given  by  C.  S.  Draper,  Notes  on  Instrument  Engineer 
ing,  M.I.T,  Instrumentation  Laboratory,  Cambridge,  Mass.,  September,  1950. 
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DERIVATION  1 

SPECIFIC  FORCE  COMPONENTS 
ASSOCIATED  WITH  INERTIAL  NAVIGATION 


The  object  of  this  derivation  is  to  obtain  the  specific  force  components 
associated  with  the  motion  of  a  point  P  representing  the  position  of  a  vehicle 
in  a  navigation  problem.  The  point  is  located,  as  shown  in  Figure  1-1,  at  a 
vector  distance  Rmp  from  the  origin  m  of  moving  coordinate  system  m. 

Origin  m  is  located  at  a  vector  distance  Rrm  from  the  origin  r  of  reference 


ORIGIN  OF 
REFERENCE 
SYSTEM 


ANGULAR  VELOCITY  OF  MOVING  SYSTEM  WITH 
RESPECT  TO  REFERENCE  SYSTEM 


Figure  1-1.  Geometrical  Relationships  of  Location  Vectors  in  Two  Coordinate  Systems. 

coordinate  system  r.  System  m  has  motion  of  both  translation  and  rotation 
with  respect  to  system  r.  A  gravitational  field  of  intensity  G  is  present  at 
point  P. 

The  vector  distance  Rfp  from  origin  r  to  point  P  is 


R  „  =  R  +  R  n 
rP  rm  mP 


Differentiation  of  eq.  (1-1)  with  respect  to  system  r  gives 


[R  „]  =  [R  ]  +  [R  J 

1  rP  r  rm  r  mPJr 


(1-D. 

d-2) 


where  f  ]  denotes  differentiation  with  respect  to  system  r,  and  (') 
denotes  d/dt.  Application  of  the  equation  of  Coriolis  to  [Rmp]r  gives 


(1-3) 
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where  W.  is  the  angular  velocity  of  system  m  with  respect  to  system  r. 
Substitution  of  this  relationship  into  eq,  (1-2)  gives 

!"rp'r  =  [5rm]r  *  (Smp]m  ♦  Wrm  x  Rmp  (1-4) 

In  words,  eq.  (1-4)  states 


velocity  of 
point  P  with 
respect  to 
system  r 


velocity  of 


velocity  of  cross  product  of  angu- 


point  m  with  point  P  with  lar  velocity  of  system 
respect  to  +  respect  to  m  with  respect  to  sys- 


system  r 


system  m  tern  rand  distance  of 
point  P  from  point  m 


Differentiation  of  eq.  (1-4)  with  respect  to  system  r  gives 


[R  _]  =  [R  ]  +  [R  n] 
1  rPJr  rnrr  mPJ 


+  W  x  R  +  W  x  [r  ] 
m  r  rm  mP  rm  mP  r 


Application  of  the  equation  of  Coriolis  to 


gives 


^RmP^m  r  =  ^RmP^m  +  Wrm  X  ^RmP^m 


Substitution  of  eqs.  (1-3)  and  (1-6)  into  eq.  (1-5)  gives 
******  * 

[R =  [Rr„]r  +  [RmDL  +  W  X  (w  XR  ]  +  W  x  R  _ 
rP  r  rm  r  mP  m  rm  rm  mP  rm  mP 

• 

+  2  W  x  [R  n]  fi¬ 
rm  mPJm  ' 


In  words,  eq.  (1-7)  states 

acceleration  acceleration 

of  point  P  of  point  m 

with  respect  -  with  respect 

to  system  r  to  system  r 

tangential  acceleration  of 

point  P  due^o  angular  ac- 

+  celeration  W  of  system 
rm 

m  with  respect  to  system  r 


,  ,.  centripetal  accelera- 

acceleration  ..  f  .  ,  _  . 
nmnf  d  tion  of  P°int  p  dui  to 


of  point  P 


+  angular  velocity  W 


with  respect  &  ^  rm 

to  System  m  of  system  m  with  re¬ 
spect  to  system  r 

Coriolis  accelera¬ 
tion  of  point  P 
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For  application  of  true  position  with  respect  to  the  Earth  let 

Point  r  in  eq.  (1-7) - ►  E,  i.e.,  center  of  Earth 

Point  m  in  eq.  (1-7) - ►  E,  i.e.,  center  of  Earth 

System  r  in  eq.  (1-7) - ►  E,  i.e.,  Earth  reference  axes 

System  m  in  eq.  (1-7) - ►  P^  i.e..  true  position  axes 


Then  eq.  (1-7)  becomes 

^REP^E  =  ^EPt^Pt  +  WEPt  X  ^EPt  X  ^EP^  +  ^EPt  X  REPt 


+  2WEPtX  ^EP^Pt 


(1-9) 


since  fREE]E  =  °- 

For  application  to  the  acceleration  of  true  position  Pt  with  respect  to 
inertial  space,  via  the  Earth,  let 


Point  r  in  eq.  (1-7)  — 
Point  m  in  eq.  (1-7)  — 
System  r  in  eq.  (1-7)- 
System  m  in  eq.  (1-7)- 


E,  i.e.,  center  of  Earth 
E,  i.e.,  center  of  Earth 
I,  i.e.,  inertial  space  reference  axes 
E.  i.e.,  Earth  reference  axes 


Then  eq.  (1-7)  becomes,  using  the  relationship  of  eq.  (1-3), 
(fiEPt  ”  [ftEP()E  +  WIE  x  (WIE  *  5EP,1  +  2WIE  x  ^EP^P, 
,2WiEx[WEptXREpt] 

since  [RggJj  =  0  and  WJE  =  0. 

Substitution  of  eq.  (1-8)  into  eq.  (1-9)  gives 

^EP^I  =  tREP^Pt  +  WEPt  X  ^WEPt  X  REpJ  +  WEPt  X  REPt 


(1-9) 


+  2W. 


EP. 


X  ^REP^Pt  +  WIE  X  ^WIE  X  REP^  +  2WIE  X  fREP^Pt 

(1-10) 


+  2WIE  X  fWEPt  X  REPtJ 
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In  marine  operation,  where  radial  acceleration  of  the  vehicle  is  small  com¬ 
pared  to  gravitational  acceleration,  and  the  radial  velocity  of  the  vehicle  is 
small  compared  to  its  surface  velocity,  [Rgp  ]p  and  [R£p  ]p  nnay  be 

neglected;  eq.  (1-10)  then  reduces  to 

[S,E  JI  =  WEP  x  lWEP  x  BEJ  +  *EP  x  5E  +  *IE  X  [*IE  X  SE 1 


+  2Wiex[wepxre1 


(1-11) 


Ep  — Rg  and  the  subscript  t  is  dropped  for  convenience,  since 


where  R 

no  ambiguity  is  possible  by  the  simplification 


Expressed  in  terms  of  the  associated  inertia  reaction  specific  force 
terms,  eq.  (1-11)  becomes 


(sf)(IR)  =  (sf)(ffi)WEp+  {Sf)U/.-Ep 


+  (sf). 


+  (sf). 


(IRjW^  T  '^'(IR)Wie  t  VkJi/(IR)C 


(1-12) 


where 


(s,,«R) 

[re]i 

(st)0R)wEp 

=  -  WEp  x  [WEp  X  Re]  =  -  WEp  x  vEp  ;  negligible 

compared  to  (sf)„  below,  at  marine  speeds, 
u 

(st)(m)wEp 

=  -^EPX«E 

(s,)(m)wIE 

=  -W  x  (W  x  R  '] ;  effect  of  Earth's  daily  rotation. 
Ik  Ik  k 

(sJ)(m)c 

=  -  2  WIE  x  [ WEp  x  Re  ]  =  -  2  WIE  X  7ep  ;  horizontal 

transverse  component  of  Coriolis  acceleration; 
negligible*  at  marine  speeds. 

*  The  magnitude  of  the  tangential  component  of  (sl)^^  is 

I^IE^EpI  (tong)  =  2vEPWIESin<Lat) 

For  Lat  =  45°  and  v^p  = 

10  knots,  the  error  in  true  vertical  indication  attributable  to  the 
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When  the  gravitational  specific  force  is  included,  eq.  (1-12)  becomes 


<st>(res)  =  (sf)G  +  (st)(IR) 


where 


(sf) 


(res) 


(sf). 


'g  '  'G 
It  is  to  be  noted  that 


=  lct\  j.  f  fgO  _  fgfi,  \  1. 

'  'S  ‘  tv  '(IRh  'bI/(IR)WIEJ 


=  resultant  specific  force 
=  G  =  gravitational  field  intensity 
=  (5f)„  +  (sf) 


(1-13) 

(1-14) 


(IR)W.„  =  g  =  gravity  field  intensity 


(sf)  is  vertical,  by  definition 
(ii)(IR)WEp  13  VerliCaI 


(sf) 


(IR)WEp  is  horizontal 


(sf)^jpj^  is  horizontal  and  transverse 


The  vector  (sf)(res)  is  the  quantity  to  which  a  pendulous  element  responds. 


tangential  Coriolis  specific  force  is 

((E)Vt)p  =  — s  "I"  m*nute  °f  arc 
This  is  the  reason  (sf)^^  is  neglected  in  the  detailed  system  discussions  to  follow. 
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MAlUCMVfim  MrttnuTl  C*  IK«K*C»0' 

tNlUtfWiM'AttON  lAtCAAlOlf 

***  f.  £  .V.  S/ 


HORIZONTAL  INERTIA  REACTION  VERTICAL  INERTIA  REACTION 
SPECIFIC  FORCE  SPECIFIC  FORCE 


Geometrical  Relationships 


s(nA(l  .  3)  .  sln|[(C)v](l>1)  .  A(, 

.  a))*  sln[(C)v](u)C0SA  „  _  aJ  .  cos[(C)V](ll)  slnA(, 

.  <s,)(in)H 

S,nV»>-  TT-i— 

(s!)<mm  *  RE(Vcrl)t 

(si)  .  <sO(1R)v 

(s(V)v°-re<v",)2. 

cosv*>* 

(Sl)g  •  B 

Note  that  the  inertia  reaction  specific  force  associated  with  the  centripetal  acceleration 
fue  to  the  Earth’s  daily  rotation  is  included  in  the  specific  force  of  gravity. 
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Figure  2-,<  Geometrical  Relationships  Associated  with  a  Pendulous  Element  as  an 

Indicator  of  the  Vertical. 


DERIVATION  2 

CHARACTERISTICS  OF  A  PENDULOUS  ELEMENT 
AS  AN  INDICATOR  OF  THE  VERTICAL 


Newton's  law  of  motion  applied  to  a  pendulous  element  that  is  tracking 
the  resultant  specific  force  is 


7  •  • 

m  r“  (Vert).  =  m  1  (sf),  >  sin  A,.  , 

p  gp  x  PP  (res)  (l-a) 


(2-1) 


where 


m  =  mass  of  pendulous  element 
P 


r^  =  radius  of  gyration  of  pendulous  element 

For  other  definitions  see  Fig.  2-1. 

Substitution  in  eq.  (2-1)  of  the  relationships  for  (Vert),  and 

(sf)(res)sxn  A^.  ^  listed  in  Fig.  2-1  gives,  with  cancellation  of 

m  and  division  by  r2  , 

P  £P 


(Vert)t  -  [(C)V]M  =  l(sf)g  +  (sf)(m)v)  sin  [(C)V](M) 

rgP 

*  Hr  “s  f<C)V)(V)  (2-2) 

r  gp 

Substitution  in  eq.  (2-2)  of  the  relationships  for  the  specific  force  com 
ponents  listed  in  Fig.  2-1,  using  the  approximations*  that  sin  [(C)VL. 

=  f (C)v aixd  cos  [(C)V]^t  j  =  1,  gives:  ’ 


*  With  any  pendulous  element  that  is  a  suitable  indicator  of  the  vertical,  the  angle  [(C)v]^  ^  should 
be  only  a  lew  minutes  of  arc,  or  at  most  a  few  degrees.  ’ 
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(Vert),  -  E(C)V]{t  i) 


V  Ee  -  re  <Vert>2t'  f<c)vl(t,D 

rgP 


+  J-  Re  (Vert), 

rgP 


(2-3) 


A  regrouping  of  the  terms  in  eq.  (2-3)  gives 


[(C)V] 


[g  -  Re  (Vert)2,]  E(C)V](y) 


1 

=  [1 - RE](Vert>t 

(2-4) 


The  condition  of  Schuler  tuning  is  that  the  coefficient  on  the  right-hand  side 
of  eq.  (2-4)  be  zero.  In  an  attempt  to  accomplish  this  let 


r  E(set) 

gP 

Substitution  of  the  relationship  of  eq.  (2-5)  into  eq.  (2-4)  gives 


[(C)V] 


(t,i)  + 


g  -  Re  (Vert)2 
RE(set) 


[<OV]<t ,i)  -  t1  - 


be 

RE(set) 


](Vert)t 

(2-6) 


DERIVATION  3 

VERTICAL  INDICATION 


A.  Basic  Principles 

The  vector  angle  representing  the  correction  to  the  indicated  vertical 
with  respect  to  the  true  vertical  is  expressed  mathematically  as 

[<^1(1,1)  =TV,  •  l3-‘> 

where  the  quantities  are  defined  in  Fig.  3-1. 

Indication  with  Respect  to  Earth  Axes  —  Accelerometer  Tracking 

The  rate  of  change  of  this  angle  with  respect  to  indicated  position  axes, 
i.e.,  in  a  set  of  coordinates  that  consists  of  the  indicated  vertical  and  two  other 
mutually  orthogonal  axes  arbitrarily  oriented  in  the  indicated  horizontal 
plane,  is 

="(EV)t  -Vlt  ->P)V1  (3-2) 

> 

where  [(C) V].,  =  rate  of  change  of  indicated  vertical  correction 

'  ’  '  i  with  respect  to  indicated  position  axes 

W/Tr,,u  =  angular  velocity  of  true  vertical  with  respect  to 
(EV)t  the  Earth 

W,„  . .  =  angular  velocity  of  indicated  vertical  with  respect 

'  '  to  the  Earth 

WfEPiV  =  anSular  velocity  of  indicated  position  axes  with 
i  respect  to  the  Earth  about  the  indicated  vertical 

The  magnitude  of  the  indicated  horizontal  component  of  the  resultant 
specific  forces,  i.  e. ,  the  quantity  sensed  by  an  accelerometer,  is  ex¬ 
pressed  mathematically  in  eq.  (3-3). 
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(res)H. 


^(res)H.  ”  H'.  *  ^S^(res) 


where 


(si),  =  indicated  horizontal  component  of  resultant 

'res’Hi  specific  force 

I  I  =  denotes  magnitude  of  the  enclosed  vector  quantity. 

Let  the  system  operate  such  that  the  angular  velocity  of  the  in¬ 
dicated  vertical  with  respect  to  the  Earth  is  a  function  of  the  magnitude 
of  the  indicated  horizontal  component  of  the  resultant  specific  force  and 
directed  so  as  to  rotate  the  indicated  vertical  toward  the  resultant 
specific  force  vector,  namely 

>V)i=<PP)(Vi)(I,W)  IV  *<sl>(res)l  <3-"> 


where 


(PF),v-wf  w>  =  vertical-indicating  performance  function  for  a 
'  ’  '  specific  force  input  and  an  angular  velocity 

output. 

Substitution  of  eq.  (3-4)  into  eq.  (3-2),  plus  the  relationships  for 
the  resultant  specific  force  from  Derivation  1,  gives 


[(C)V](t,»P.  =W(EV)t  -‘PfW,W)TV,  X{ 


g  “  ^WEP^E  X  RE 


"  WEP  X  <WEP  X  V  -  2WIE  X  ^WEP  X  M  "  W(EP)V.  X  t(C)V^(t,  i) 


The  specific  force  components  in  eq.  (3-5)  can  be  conveniently 
resolved  into  components  parallel  to  the  true  vertical  and  perpendicular 
to  the  true  vertical  (the  arbitrary  nature  of  axes  in  the  horizontal  plane 
requires  no  further  subdivision  into  components),  namely 

i  =iv  g 

t 


RE  =  !vt  RE 


WEP  "  \  W(EP)Vt  +  W(EV)t 


MASSACHUSETIS  INSIIIUU  Of  UCHN010G* 
INSIRUMENlATjON  OftORATOXf 
app.  w  ££%Ht 


VECTOR  FORM 

OF  [(C)V](M) 
[<C)V](lii) 


RESULTANT 
SPECIFIC  FORCE 


UNIT  VECTOR  ALONG 
INDICATED  VERTICAL 


UNIT  VECTOR  ALONG 
TRUE  VERTICAL 


'(EV)t  \  / 
ANGULAR  VELOCITY 
OF  TRUE  VERTICAL 
WITH  RESPECT 
TO  EARTH  j 


APPARENT 

VERTICAL 


INDICATED 

VERTICAL 


TRUE' 
VERTICAL 


(t,L) 

CORRECTION  TO  INDICATED 
VERTICAL  WITH  RESPECT  TO 
TRUE  VERTICAL.  (Angle  through 
which  Indicated  Vertical  must 
be  rotated  to  align  it  with 
True  Vertical.) 


AS  A  VECTOR  ANGLE 

[(£)v]<‘.»  *  \  x  'Sc 


(U)  s  7  x  T 


t,i 


,  SINCE  [(C)v](t>i)  IS  SMALL 


Figure  3-1.  Geometrical  Relationships  Associated  with  the  Vector  Angle  Form  of 
the  Correction  to  the  Indicated  Vertical  with  Respect  to  the  True  Vertical. 


62 


Substitution  of  these  relationships  into  eq.  (3-5)  gives  the  specific  force  in 
terms  of  components  along  the  true  vertical,  along  the  angular  velocity 
vector  of  the  true  vertical  (which  must  be  perpendicular  to  the  true  vertical) 
and  perpendicular  to  both  of  these,  namely 

KC*Vl(t,  i)P.  =  W(EV)t  “  !PF)(Vi)(f,  W)V  X  {  V,  ^  ~  REW2(EV)t 

-  2V*ie>  •  2(Wie  •  tyi 

+  [(«EP>E  x  \]Re}  -  V>>V,  X  KW](l>  „  (3-6) 

Carrying  out  the  vector  operation  indicated  in  eq.  (3-6),  regrouping 
terms  and  omitting  the  components  parallel  to  the  indicated  vertical  as 
not  being  applicable  to  the  indication  of  the  vertical,  gives 

KC)Vl(t,i)P.  +  (PF)(Vi)(f,W)  (sf)(res)Vt  [(c>v](t,  i)  =  W(EV)t 

-*PFVi)(f,W)IVW(EV)t*P1  “  *PF*(Vi)(t,W)RElW(EP)Vt 

+2(SW  \»  ]  <TV.  x  W(EV>t>  -  W(EP)V.  x  D  (3-,) 

where 

^sR{res)V,  -  6  *  REW2(EV)t  '  2RE»IE  '  W(EV)t>  <3-8> 

=  component  of  resultant  specific  force  along  the  true 
vertical 

Let  (PF)^..^  be  the  following  operator: 

(PF)(Vi)(f,W)  =  S(Vi)(f,W)  Jp  (  )dt  (3"9) 

i 

where  /_  denotes  integration  of  a  function  of  the  indicated  axes  with  re- 
i 

spect  to  the  time.  Substitution  of  equation  (3-9)  into  eq.  (3-7)  gives  eq.  (3-10): 
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'(C>Vl(t,i)Pi  *  S(Vi)(f,W)  f  <SI)(res)Vt  f<C1V)(t,i)  *  '  W(EV)t 
1  i 

"  S(Vi)(f,W)  fp  RE(^(EV)t)P.dt  '  S(Vi)(f,W)  fp  RE[W(EP)Vt 
i  i 

+  2(w1E-TVt)](Tv  X  W(EV),,dl  -  w  X  imM  (3-10) 

Differentiating  (3-10)  with  respect  to  the  time  gives  the  differential  per¬ 
formance  equation  in  terms  of  functions  of  the  indicated  axes: 

=  d  -  S(Vi)(f,W)RE)(W(EV)t)P.  ■  S(Vi)(f,W)RE^(EP)Vt 

C - * - - v 

+  2  TVt)]  (1V.  X  W(EV)t>  -  ^(EPjV.  X  t<C>V  WPi  {3_11) 

where  it  is  assumed  that  variations  in  Rg  are  slow  enough  to  be  treated 
as  quasi-static  and  variations  in  (sf)(res)y  are  either  slow  enough  to  be 

quasi-static  or  fast  enough  to  be  ineffective. 

For  Schuler  tuning  let 

S  ♦  -  -!_ 

tVi)(f,W)  “  Rs 

where  R  is  a  set  value  made  as  nearly  equal  as  possible  to  the  magnitude 
of  R£.  If  it  is  then  assumed  that  Rg  =  Rg,  eq.  (3-11)  becomes 

.  ••  ^(res)V.  _  _  —  _  _ 

[(c)vJ(t  i)pi  +  C(c)v](t  i)  =  -  fW(EP)Vt  +  2(W(IE)*1Vt^(1V. 


X  W(EV)t>  -  ^(EPJV.  X  I<C)VWP| 


(3-12) 


The  terms  W^Ep^  (ly  x  W(ev)P’  and  2  ^Ie’  1V^^Vi  X  W(EV)t^  are 


(res)V. 


small  in  comparison  with 
shown  as  follows: 


,  und  can  be  neglected.  This  can  be 
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W(EP)V  ’  m:>onirade  °*  the  ra*e  °f  change  of  the  meridian  about  the 
true  vertical  with  re&pnut  to  the  Earth,  is  given  by 


W(FP)V  =  IT-  sin(c°)tan  (La1) 
t  E 


(3-13) 


where 


Vg  -  Velocity  of  base 
Co  =  course  of  base 
Lat  =  latitude 


(EV)t  Rt 


(3-14) 


r(EP)Vt  (1V.  X  W(EV)tT  W(EP)VtW(EV)t 


(3-15) 


2  (WTP.  lv  )(lw  x  W/1?„j|  3  2  W(Ey)t  WIE  sin  (Lat) 


IE  (EV)t 


Assume  that  the  base  is  traveling  due  east  at  10  knots  at  70°  latitude: 


w  w 
(EP)Vt  (EV)t 


;(res)V. 


RsW(EP)VtW(EV)t 


(res)V. 


=  00.02  min.  of  arc 


(3-16) 


2Wlg  sin  (Lat)  W(Ey)t  2RgWIE  sin  (Lat)  W(Ey)t 


'(res)V. 


;(res)V, 


00.25  min  of  arc 
(3-17) 


Neglecting  the  above  terms,  eq.  (3-12)  reduces  to 


[(c)vl(t,i)P. + 


'(res)V, 


«C>V](M)  ’  -  (W(EP)V.  X  [<C>Vl(t,t)> 

(3-18) 


The  term  O^gp^y  X  [(C)V]^.  .^)  represents  an  orientation  transfer,  not 
a  magnitude  effect,  and  represents  a  secondary  coupling  between  the  x  and 
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y  systems.  Since  this  will  only  produce  a  "spiraling  in"  effect  as  the 
system  approaches  its  null  position,  it  will  be  ignored  in  subsequent 
discussions. 


Equation  (3-18)  now  reduces  to 

(s"(res)  V,  - 

[(C>Vi(t,iV— - ^(M)  *° 

i  s 


(3-19) 


B.  Instrumentation  of  a  Single-Axis  System 
Definitions  of  Symbols  Used  in  the  Following  Derivations 


EB 


IB 


P. 


(Pi>o 


[(c)P](t,i) 

^P](t,i)o 


f(C)V] 

[(C)V] 

[(C)V] 


'( ) 


[(C)V] 
e 
G 

g 

H 

l(  ) 


(t,i) 
(t,  i)o 
(t,r) 
(t,  r)o 


Coriolis  acceleration 

Acceleration  of  base  with  respect  to  Earth 
Acceleration  of  base  with  respect  to  inertial  space 
Indicated  position  (latitude  or  longitude)  of  vehicle 
Indicated  position  of  vehicle  at  start  of  problem 
True  position  of  base 
Position  correction 

Position  correction  at  start  of  problem 

Correction  to  the  indicated  verti,  1 

Correction  to  indicated  vertical  at  start  of  problem 

Correction  to  reference  vertical 

Correction  to  reference  vertical  at  start  of  problem 

Output  signal  voltage  of  device  indicated  by  subscript 

Gravitational  field  vector 

Gravity  field  vector 

Angular  momentum  of  gyro  unit  rotor 

Output  current  of  device  indicated  by  subscript 
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M 

g 

Torque  on  gyro  gimbal 

r  =  re 

Radius  of  the  Earth,  a  function  of  geographic  position 

R 

s 

Set  value  of  Earth  radius  used  for  Schuler  tuning 

S(  )(  ,  ) 

Sensitivity  of  component  or  component-group  indicated 
by  first  subscript,  with  input  and  output  indicated  by 
second  subscript 

V. 

l 

Indicated  vertical 

(V.) 

10 

Indicated  vertical  at  start  of  problem 

V 

r 

Reference  vertical 

Vt 

True  vertical 

(Vo 

True  vertical  at  start  of  problem 

W(EV)i 

Angular  velocity  of  indicated  vertical  with  respect  to 
the  Earth 

W(EV)r 

Angular  velocity  of  reference  vertical  with  respect  to 
the  Earth 

W(EV)t 

Angular  velocity  of  true  vertical  with  respect  to  the 
Earth 

WIB 

Angular  velocity  of  base  with  respect  to  inertial  space 

WIE 

Angular  velocity  of  Earth  with  respect  to  inertial  space 

W(IV)i 

Angular  velocity  of  indicated  vertical  with  respect  to 
inertial  space 

cmds 

Controlled  member  drive  system 

Figure  3-2  is  a  mathematical  functional  diagram  of  a  single -axis  system 
designed  to  indicate  the  vertical  under  dynamic  conditions.  The  component 
will  be  considered  as  ideal  and  the  system  to  be  undamped  (non-ideal  com¬ 
ponents  and  damped  systems  will  be  analyzed  in  the  second  report),  and  it 
will  be  shown  that  the  equation  of  motion  of  the  system  is  (3-19). 

The  quantity  [(C)Vj^.  ^  is  defined  by  equation  (3-1);  that  is,  [(C)V]^t 
is  the  angle  between  the  true  vertical  and  the  indicated  vertical. 
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Then,  in  scalar  form  and  operator  notation 

p[<c)V](M)  -  w(EV)t  -  w(EV)i 

where  p  =  ^- 

Equation  (3-20)  corresponds  to  eq.  (3-2)  with  the  coupling  term 

<Vp)vi  neglected- 

The  proposed  system  consists  essentially  of  an  accelerometer,  an 
integrator,  a  single-degree-of -freedom  integrating  gyro  unit,  and  a 
drive  system  for  the  stable  platform  or  controlled  member.  The 
accelerometer  is  mounted  on,  and  near  the  axis  of  rotation  of,  the 
controlled  member.  It  will  therefore  sense  the  indicated  acceleration 
of  the  base  with  respect  to  inertial  space,  or  the  projection  of  the 
resultant  specific  inertia  reaction  force  on  the  indicated  horizontal 
plane  (in  general  this  force  will  contain  both  acceleration  and  gravity 
effects). 

The  performance  equation  for  the  accelerometer  unit  is  (neglecting 
accelerometer  dynamics*): 

e(au)  =  S(au)(a,  e)^C)V^(t,  i)  +  ^B***  (3’22) 

where  e&u  =  accelerometer  unit  output  signal 

g 

(au)(a,  e)  =  accelerometer  unit  sensitivity  for  an  acceleration 
input  and  a  voltage  output 

a.„  =  acceleration  of  the  base  with  respect  to  inertial  space 

Id 

Since  the  system  intrinsically  minimizes  [(C)VL.  .,,this  angle  will  be 
small  enough  in  general  to  permit  the  assumptions:  '  ’ 

sm  [(C)v](t  t)  =  l(C)v](t  „ 

c°s[(C)V]((l)3  1 

In  Derivation  1  it  was  shown  that,  at  submarine  velocities 

alB  "  R»W(EV)t  (3-23) 


*  To  incorporate  Schuler  tuning  characteristics,  the  system  is  designed  to  have  a  natural  period 
ot  approximately  84  minutes.  Therefore,  the  dynamics  associated  with  the  accelerometer  unit 
and  the  controlled  member  drive  are  neglected. 

**  In  this  scalar  equation  the  centripetal  acceleration  due  to  the  Earth's  daily  rotation  is  included 
in  g  rather  than  ajQ  .  The  net  vertical  force  is  effectively  g  for  submarine  operation. 


(3-20) 

(3-21) 
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Substituting  this  in  (3-22)  gives 


i(au)  =  S(au)(a,e)W«V](M)+KpW(EV)t) 


(3-24) 


The  integrator  performance  equation  is  given  by 
*(int)  “  S(int)(e,i)  "p"  e(au)  +  *(int)o 


(3-25) 


i(int)  =  current  outPut  integrator 


5(int)(e,i)  =  sensitivity  of  integrator  for  voltage  input,  current  rate  output 


(int)o  =  integrator  current  at  start  of  problem 

The  primary  purpose  of  the  gyro  unit  is  to  provide  a  means  of  stabi¬ 
lizing  the  controlled  member  against  arbitrary  motions  of  the  base;  however, 
in  the  present  analysis,  the  gyro  unit  is  simply  considered  to  be  part  of  the 
angular  velocity  generating  system.  The  input  current  produces  a  torque  on 
the  gyro  unit  which  is  balanced  by  the  angular  velocity  of  the  controlled  member 
with  respect  to  inertial. space,  thereby  keeping  the  gyro  unit  on  null, 

S(tg)(i,  M)^(int)  +i-]  =  HW—  =  HW—  <3-26> 


iw]=HW 


(cmds)  =  HW(IV)i 


where 


=  sensitivity  of  gyro  unit  torque  generator  for  current 


(S)(h  I  inpUt  and  torque  output 

i  =  current  compensating  for  rotation  of  earth  with  respect  to 

inertial  space 

H  =  angular  momentum  of  gyro  unit 

W 

(cmds)  =  angular  velocity  of  controlled  member  with  respect  to  inertial 
space 

W  £  w 

(IV)i  ~  (cmds)  =  angular  velocity  of  indicated  vertical  with  respect  to 
inertial  space 


5(cmds)(i,W) 


_  (tg)(i,M) 


(3-27) 


where 


3(cmds)(i,  W) 


sensitivity  of  controlled  member  drive 
system  for  a  current  input  and  an  angu¬ 
lar  velocity  output 


From  the  foregoing, 


It 

> 

W(EV)i+  WIE  "  S(cmds)(i,W)i(int)  +  S(cmds)(i,  W)V 

(3-28) 

Set 

.  .  w* 

S(cmds)(i,  W) 

(3-29) 

then 

W(EV)i  =  S(cmds)(i,W)1(int) 

(3-30) 

W(EV)i  “  S(au)(a,  e)  S(int)(e,  i)S(cmds)(i,W)  p  (g[(C)V](t,  +  RpW(Ey)t) 

+  S(cmds)(i,W)i(int)o  (3-31) 

Equation  (3-31)  corresponds  to  equation  (3-4). 

Substituting  eqs.  (3-31)  into  (3-20)  and  differentiating  the  result, 

[p2  +  S(au)(a, e)  S(int)(e,i)S(cmds)(i,W) ^C)V^(t,i) 

“  ^  ”  S(au)(a,  e)S(int)(e,  i)  S(cmds)(i,W)R^  P  W  (EV)t  (3-32) 
For  Schuler  tuning,  set 

S(au)(a,e)S(int)(e,i)S(cmds)(i,W)=:  R~  (3'33) 

s 

Eq.  (3-32)  reduces  to 

(p2+  ~K~  )^(C)Vl(t,i)  =  °  (3'34) 

s 

and  is  identical  with  (3-19). 

The  solution  to  (3-34)  is 

■  KC)V)(t,i)(<S)  oosySIt.  t<C)V1(M)0^r  sinyEt 

(3-35) 

where 

[(C)V](t  =  Correction  to  indicated  vertical  at  start  of  problem 

t(C)VL  ..  =  Rate  of  change  of  correction  to  indicated  vertical  at 
’  start  of  problem. 
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gyro  unit,  i.e.,the  ratio  of  the  sensitivity 
lingular  velocity  input-  torque  output,  these 
It  controlled  member  drive  is  suf  ficiently  fast 
main  feedback  loop. 


’(cmds)  (L,W) 
P 


(cmds) 

('SR), 


‘1 


(gu)UtW) 


Vi 


I 


(gu) 

ccgu) 

CONTROLLED 

MEMBER 

SYSTEM 

AND  DRIVE 

* 

Vi 


■■  Meevy  fine  indicate s  rigid 
mechanical  connections 

—  Light  fine  indicates  signal  transfer 

The  cent  rife  get  force  due  to  the  earth's 
dally  rotation  affects  the  accelerometer 
through  the  support.  (Set  Fig.  l.J 


MAUACHUIITTS  INSTITUTE  OF  TtCHNOiOOT 
INSTRUMMTATION  UIOMTOIY 

*»■  *  **  »•  /*•  x/ 


Figure  3-2.  Mathematical  Functional  Diagram  of  a  Basic  Single-Axis  Vertical  In¬ 
dicating  System  Using  the  Components  Shown  in  Figure  1 5. 


DERIVATION  4 


POSITION  INDICATION  BY  MEANS  OF 
OPEN-CHAIN  INTEGRATION  OF 
THE  ANGULAR  VELOCITY  OF  THE  INDICATED  VERTICAL 


Figure  4-1  is  a  mathematical  functional  diagram  of  a  single-axis  system 
for  the  indication  of  position.  Here,  the  input  to  the  gyro  unit  in  Derivation  3B 
is  integrated  once  to  give  the  indicated  position.  The  gyro  unit  input  is  pro¬ 
portional  to  the  angular  velocity  of  the  indicated  vertical.  Thus,  the  indicated 
position  is  given  by: 


Pi  =  S(Pi)(i,P)  p  l(int)  +  (Pi)o 


(4-1) 


where  (P.)  is  the  indicated  position  at  the  start  of  the  problem. 

From  (3-30),  i^.^  is  related  to  the  angular  velocity  of  the  indicated  ver¬ 
tical  as  follows: 


W 


(int) 


(EV)i 


(cmds)(i,W) 


[3-30] 


and  from  (3-21) 


P«C)V](t ,i)  =  W(EV)t  "  W(EV)i 


[3-20] 


W 


P,  =  S 


(EV)t 


-  S, 


1  p[{C)V](ti) 

m  4*  5 - M  +  (P. 


i  °(Pi)(i,P)  P  S(cmds)(i>W)  (Pi)(i,P)  p  S(cmds)(ijW)  '  ^(4-2) 


Set  the  magnitude  of  the  position  integrator  sensitivity  equal  to  that  of  the 
controlled  member  drive  system  sensitivity: 

S(Pi)(i,P)  =  S(cmds)(i,W)  (4_3) 

(This  constitutes  an  arbitrary  adjustment  of  the  system  parameters) 

Note  that 
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Figure  4-1.  Mathematical  Functional  Diagram  of  a  System  for  Indicating  Position 
by  Means  of  an  Open-Chain  Integration  of  the  Angular  Velocity  of  the  Indicated 
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DERIVATION  5 


POSITION  INDICATION  BY  MEANS  OF 
OPEN-CHAIN  INTEGRATION  OF  THE 
ANGULAR  ACCELERATION  OF  THE  INDICATED  VERTICAL 


Figure  5-1  is  a  mathematical  functional  diagram  of  another  system  for 
position  indication.  Here,  the  accelerometer  unit  output  in  the  basic 
vertical  indicator  is  integrated  twice  with  respect  to  the  time  to  give  the  in¬ 
dicated  position.  The  accelerometer  unit  output  is  proportional  to  the  angu¬ 
lar  acceleration  of  the  indicated  vertical;  this  is  discussed  in  the  foregoing 
text.  The  indicated  position,  from  Fig.  5-1,  is  given  by: 

P.  =  S/n  n\  “V  ei  \  +  (P-)  t  +  (P.)  (5-1) 

l  (Pdi)(e,P)  ^2  (au)  l  o  i  o 

From  (3-22)  and  (3-23)  in  Derivation  3  the  accelerometer  unit  output  is: 

e(au)  =  S(au)(a,e)fe!<C)Vl(t,i)  +  RpW(EV)t>  (5-2> 
Adjust  the  system  parameters  so  that 


S(Pdi)(e,P)  -  RsS(au)(ae) 

(5-3) 

p  -  —  & 

l  2  R 

p  s 

From  (3-34) 

[<c>v](M,  ♦  -V  Pw(EV)t  ♦  (P.)ot  +  (P.)o 

(5-4) 

<P2  +  T-  >  "C)Vl(t,i)  *  0 

[3-34] 

J  “2  f<c>vI<t,i)  -  t<c>v](M>  -  t(C)Vl(U)o  *  t<C’)VW  (5-5) 
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pi  *  -  I<c>vl(t,i> >  f<c>vJ(t,i)o  *  f<c>v)(M)ot  +  pt  -  <pt>o 


(P.)  t  +  (P.)  t  *  (p.) 

'  t'o  10  10 


(5-7) 


The  correction  to  the  indicated  position,  from  (4-7)  in  Derivation  4,  is 

[(c)p!(M)  -  pt  -  pt  l4-1' 

Also 

[<C,Pl(t,i)o  ■  (pt>0  -  (Pi>0  [4-8) 

t<C>P](t,i)o  ■  <Pt>o  '  <Pi>o  (5-8> 

From  eq.  (3-35)  in  Derivation  3 

[(C)V]M  .  t<c>v]<t,i)o  $ 

[3-35] 

Substituting  the  above  into  (4-7)  gives 

[<c>pI(t,l)  ■  t<°>p](t,i)o  *  t<C>p](M)o‘  '  '(c)v)(t,t)0  (1  -  C°S^F‘> 


R 


^C)Vl(t,i)o  ^  sinV^~  ^ 


g  YR 


(5-9) 


Figure  5-1.  Mathematical  Functional  Diagram  of  a  System  for  Indicating  Position 
by  Means  of  an  Open-Chain  Double  Integration  of  the  Angular  Acceleration  of  the 

Indicated  Vertical. 
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DERIVATION  6 

POSITION  INDICATION  BY  MEANS  OF 
DIRECT  DOUBLE  INTEGRATION  OF  ACCELERATION 


Figure  6-1  is  a  mathematical  functional  diagram  of  a  system  for  indi¬ 
cating  position  by  the  use  of  an  accelerometer  unit  which  is  itself  doubly- 
integrating  with  respect  to  time.  Its  output  signal  is  therefore  proportional 
to  indicated  position  P^.  The  indicated  position  is,  from  Fig.  6-1: 

Pi  ~  S(ind){e,P)e(diau)  ^  ^ 

The  accelerometer  unit  output  is,  taking  into  account  the  double  integration 
action  of  the  accelerometer  unit, 

e(diau)  “  S(diau)(a,e)  2  fe[(C)V]^  .)  +  RpW(EV)P  +  ^(diaujo*  +  e(diau)o 

P  (6-2) 

To  perform  the  indicated  double  integration  in  (5-4),  it  is  necessary  to 
evaluate  [(C)  V]^  ^  for  the  vertical  indicator  used  here,  which  is  slightly 

different  from  that  of  Fig.  5-1,  since  it  involves  the  doubly -integrating 
accelerometer  unit.  This  unit  is  followed  by  a  differentiator,  whose  output 
i  iff ,  proportional  to  the  angular  velocity  of  the  indicated  vertical,  feeds  the 

gyro  unit: 

*diff  “  Sdiff(e,i)pe(diau) 

If  the  gyro  unit  is  properly  compensated  for  Earth  rate, 

^(EV)i  “  S(cmds)(i,W)'diff 
From  (3-21)  in  Derivation  3, 


P [  (c)v  ](t;i)  =  W(EV)t  "  W(EV)i 


[3-20] 
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Then 


+  RpW(EV)t)  +  e(diau)o 


p[(c)v)(t>i)  =  w(EV)t  ~  S(diau)(a,e)Sdiff(e,i)S(cmds)(i(W)  p  (gf(C)V^(t,i) 

(6-5) 

The  time  derivative  of  (6-5),  collecting  terms  in  [(C)V]^  is 

(p2  +  S(diau)(a,e)Sdiff(e,i)S(cmds)(i,W)g)  f(C)Vj(t,i) 

=  (1  '  S(diau)(a,e)Sdiff(e,i)S(cmds)(i,W)R)pW(EV)t 

(6-6) 

To  minimize  the  right  hand  side  of  (6-6),  i.e.,  to  apply  Schuler  tuning, 


let 


S(diau)(a,e)Sdiff(e,i)S(cmds)(i,W)  Rg  ^ 

and  assume  Rg  =  R  in  (6-6).  The  result  is  (3 -34), which  now  implicitly  con¬ 
tains  (6-7)  instead  of  (3-33)  as  the  Schuler  tuning  condition. 

With  this  implicit  modification,  [(C)V]^  is  given  by 

[<C)Vl(M,  -  («CV](y)o.oS^t  +  sin^t 


which  is  equivalent  to  the  result  (3-35)  of  Derivation  3. 
From  eqs.  (6-1)  and  (6-2) 

Pi  =  S(ind)(e,P)S(diau)(a,e)  ~2  (gf(C)V](t,i)  +  RpW(EV)t} 


+  S/f 


(ind)(e,P)e(diau)ot  +  S(ind)(e,P)  e(diau)o 


Now  adjust  the  system  parameters  so  that 

1 


S„ 


(ind)(e,P)  "  S(diau)(a>g) 


•v  R 


(6-8) 


(6-9) 


(6-10) 


2 


"  S(ind)(e,P)®(diau)o 


(6-11) 


<Pi>o  “  S(ind)(e,P)e(diau)o  ^  ^ 

Equation  (6-9)  becomes 

Pl4-f  +  "T  pW(EV)t  +  <Pi>o‘  +  <Pi>o  (6-13> 

p  s  '  ’  p  ' 

This  is  identical  to  equation  (5-4)  of  Derivation  5.  By  carrying  out  the 
subsequent  operations  as  given  in  that  derivation,  the  correction  to  the  indi¬ 
cated  position  is 


S(diau)(a,e) 


®(diou) 


I  (cmds) 

|  S(cmds)(i,W)  CSR)(gU)ti>W) 

!  P  =  P 


1 


CONTROLLED 

MEMBER 

•(gu) 

POWER  GENERATING 

SYSTEM 

AND  DRIVE 

u 
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Figure  6-1.  Mathematical  Functional  Diagram  of  a  System  for  Indicating  Position 
by  Means  of  a  Doubly-Integrating  Accelerometer  Unit. 


DERIVATION  7 


AN  AZIMUTH  STABILIZATION  SYSTEM 


Figure  7-1  shows  a  mathematical  functional  diagram  of  a  y-axis  vertical 
indicator  (similar  to  Fig.  3-2,  Derivation  3)  and  an  azimuth  stabilization  loop. 
Let  X,  Y,  Z,  be  the  true  position  coordinates  of  the  system  and  let  x,  y,  z,  be 
the  indicated  position  coordinates..  Coupling  between  the  y  and  z  vertical  in¬ 
dicating  loops,  ignored  in  Derivation  3,  is  a  first-order  effect  here.  Note  that 
the  Earth-rate  correction  is  not  required  in  the  y-axis  vertical  indicator. 

From  Fig.  7-1  and  Derivation  3 

e{au)y  *  S<au><a,e)MC>Vj(.,i)y  +  R»W[(EV)t]y] 


i(int)y  "  S(int)(e,i)y  p  e(au)y  +  1[(int)y]o 


(7-2) 


The  torque  balance  equation  for  the  y-gyro  unit  is 

S(tg)(i,M)yi(int)y  “  HW[(IV)i]y 
Let 

c  -  S(*g)(i>M) 

(cmds)(i,W)y  “  H 


(7-3) 


(7-4) 


The  component  of  the  angular  velocity  of  the  indicated  vertical  with  re¬ 
spect  to  inertial  space  sensed  by  the  y-gyro  unit  is  (neglecting  coupling  with 
the  x-axis  vertical  indicating  system) 

W[(IV)i]y  =  [WIE  *  (L“n)I  005  (Ut)  sin  l(C)Nl(t,i>  +  W[(EV)iJy 

(7-5) 

where  [(C)N](t  ^  is  defined  as  the  correction  to  indicated  north,  i.e.,  the 
angle  between  indicated  north  and  true  north;  and 


(Lat)  =  latitude 
(Lon)  =  longitude  rate 
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Wr i  =  y-component  of  the  angular  velocity  of  the  indicated 
*■'  ^  vertical  with  respect  to  the  Earth. 

Since  the  function  of  the  azimuth  stabilization  loop  is  to  minimize  [(C)nL,  .> 
it  is  allowable  to  set  ' >l 

sin  [(C)N](m)  S  [(C)N](t)i)  (7-6) 

From  the  definition  of  [(C)VL  .*  (Derivation  3B),  the  y-component  of  its 
time  derivative  is  '  ’  ' 

p[(C)Vl(M)y=Wt(EV,.]y-Wf(KV)i]y  <’-7> 

Substituting  (7-1), (7 -2),  (7-3),  (7-4),  (7-5)  into  (7-7),  and  differentiating  the 
result,  gives 

[p  +  S(au)(a,e)yS(int)(e,i)yS(cmds)(i,W)yg]  ^C^(t,i)y 

t1  "  S(au)(a,e)yS(int)(e,!)yS(cmds)(i,W)yR]  pW[(EV)t]y 

+  p[[WIE  +  (Ltin)]  f(C)N](t  i)  cos  (Lat)]  (7-8) 

The  Schuler  tuning  condition  is 

S(au)(a,e)yS(int).(e,i)yS(cmds)(i,W)y  ~  ~1“  ^7'9) 

o 

Equation  (7-8)  now  reduces  to 

(P2  *  -f-)  «C)V](y)y  =  [WIE  ♦  (Lin)J  t(C)N](w)  cos  (Lat) 
s 

which  is  (3-34)  except  for  the  non-zero  right-hand  side  created  by  the 
aforementioned  coupling  between  y-axis  and  z-axis  systems. 

From  Fig.  (7-1),  with  icp  =  VS(OT*)(i>w)yBs. 

‘(mult)z  -  Knt)<t,i)z  J  [l(tnt)y  +  s(cmds)(i,w)yRs  ! 

-  (t/int)x  +  iw)  sin  (Lat)j  sec  (Lat) 


(7-11) 
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where  is  the  northward  speed  of  the  ship. 

From  the  torque  balance  equation  at  the  y-gyro  unit  (see  equations  (7-3), 

(7-4)  and  (7-5)) 

(C)Nl(W)  [WIE  *  (L“n)l  cos  (Lat>  +  W[(EV)i]y] 

(7-12) 


V  ~  \int)y 


1 


(cmds)(i,W)y 


[t 


From  the  torque  balance  equation  at  the  x-gyro  unit  (neglecting  coupling  with 
the  y  and  z  systems) 


*x  ~  ^(int)x  +  *W  -  S 


(cmds)(i,W)x 


[WIEoos  (Lat)  *  W((EV)lJx] 


(7-13) 


1  =  1 

z 


S  * 

(mt);i,i)z - 1_  f  [(C)NJ  [w  +  (Lin)]  cos  (Lat) 

(mult)l  S(cmds)(l,W)y  P  1  <‘-1)  IE 


*  W[(EV,.lr^}  -<“>  '  [W^C°SM 

+  W[(EV)i]x1  «■»*>  <’-14> 

From  the  torque  balance  equation  at  the  z-gyro  unit  (neglecting  coupling  with 
the  x-axis  system) 


*z  =  S 


(cmds)(i,W)z 


[-  WIE  sln  <Lat)  +  W(EN)t  -  l(C)Vl(t,i)yWIE  cos  <Ut>] 

(7-15) 


or 

W 


Set 


-  S,  .  w.  i  +  WT„  sin  (Lat)  +  [(C)NL  W  cos  (Lat) 
(EN)i  (cmds)(i,W)z  z  IE  '  lv  (t,i)y  IE 

(7-16) 

S(cmds)(i,W)x  ”  ^(cmds)(i,W)y  ”  ®(cmds)(i,W)z  (7-l7) 

The  time  rate  of  change  of  [(C)Nj^  ^  with  respect  to  the  base  is  given  by 

»[<5m](W)  «  w(EN)l  -  w(EN)i  -  wEB  X  [<W](M) 


(7-18) 


where  W£g  =  angular  velocity  of  the  base  with  respect  to  the  Earth.  At  marine 
speeds  the  cross-product  term  may  be  neglected*.  Equation  (7-18)  now  becomes 


p[(c>N](tji)  =  W(EN)t  -  W(EN). 


(7-19) 


Substituting  equations  (7-14),  (7-16),  and  (7-17)  into  (7-19) 


p[(C)N)(M)  =  »m  -  S(int)(l  i)z  |  [  [(C)N)(m)  [Wie  +  (Lin)]  cos  (Lat) 

y 

+  W[(EV)i]y  +  TT  ]  S6C  (Lat)  +  [WIE  C0S  (Lat)  +  W[(EV)i]x]  (Lat) 


-  WIE  sin  (Lat)  -  [<C)V](t  i}  WI£  cos  (Lat) 
From  geometrical  considerations 


W[(EV)t]x  =  (L°a)  cos  <Lat) 


W(EN)t  =  “  ^Lon)  sin  (***) 


W(EN)t  "  "  W[(EV)t  Jx  ta"  (Lat) 


(7-20) 


(7-21) 


[(EV)tjy  R 


Also,  assume  that 


W[(EV)i]x  =  W[(EV)t]x 


(7-22) 


W[(EV)i]y  =  W[(EV)t3y-P[(C)V](t,i)y  (7"23) 


(7-24) 


Substituting  the  foregoing  into  (7-20)  and  differentiating, 

*  Assume  the  maximum  velocity  of  the  ship  is  50  knots.  The  order  of  magnitude  of  W„_.  is  then 

Cjo 

WEB  84.4  ft/sec  . . -6  rad 

— —  = - 2—  r~  =  4  X  10  — - 

R  20.9  X  10®  ft  SCC 

Assuming  a  correction  angle  of  10  minutes  of  arc  and  the  angle  between  W^g  and  |(C)N]^  ^ 

equal  to  90  degrees,  the  magnitude  of  the  cross  product  term  will  be  of  the  order  of  10’®. 

This  term,  then,  may  be  neglected  when  compared  with  W^ENjt  and  W^E{^(,  each  of  which 

has  a  magnitude  of  the  order  of  10’®. 


[p  +  S(int)(i,i)z  ^WIE  +  (L°n)3]  f(C)N^(t,i)  =[S(int)(i,l)z  sec  (Lat3p^(C)V^(t,i)y 
-  p{[(OV](M)y  [WIE  cos  (Lat)]j  (7-25) 

The  final  term  in  the  above  equation  is  small  compared  to  the  preceding  term 
and  can  be  neglected: 

fr2  +  S(int)(i,i)z  fWIE  +  [<C>N!(t,i)  =  [S(tnt)(i,l)zs‘i':<La‘)]p[(C)V](t>i)y 

(7-26) 

Since  the  time  variations  in  (Lon)  and  (Lat)  are  small  when  compared 
with  the  variations  in  [(C)Nj^  ^  and  [(C)V]^  ,  it  will  be  assumed  that 

(Lon)  and  (Lat)  are  quasi-static.  This  simplifies  the  subsequent  analysis. 

Then  (7-10)  and  (7-26)  combine  to  give  two  simultaneous  operator  equations 
of  the  form 


[<p4  ♦  -jE-p2  +  s(int)(1;i)z [wIE  ♦  (Lin)]  J-]  [(C)V](t  l)y  =  0 
s  s 

[(p4  +  -jf-p2  *  S(lnt)(l  i)z  [WIE  ♦  (Lin)]  -§-]  t(C)N](  J  =  0 
s  s 

The  most  general  solution  for  either  [(C)V]^  ^y  or  [(C)Nj^  ^  is 

r  ,  Pi1  p2 4  P3 1  P4t 

[(C)N](t,i)  =  Cle  +°2e  +  V  +C4e 


(7-27) 

(7-28) 


(7-29) 


where  Cj,  C g,  Cg  and  C4  are  constants  and  p^,  Pg,  Pg  and  p^  are  the  var¬ 
ious  solutions  of 


4 

P 


+  ^(int)(i,i)z^WIE 


+  (Lon)]  -~- 
s 


=  0 


(7-30) 


That  is, 


4S(int)(i,i)z^WIE  +  ^n)3Rs 
g 


(7-31) 


where  j 


The  system  will  be  stable  (i.e.,  oscillatory  with  limited  amplitude)  only  if  the 
quantity  under  the  inner  radical  is  positive;  i.e.,  if 


4S 


(int)(i,i)z 1  IE 


[WT_  +  (Lon)] 


g 


<  1 


(7-32) 


or 


S(int)(i,i)z 


4[WIEg+  (Lon)  J  R£ 


(7-33) 


Assume  that  this  condition  is  met  by  suitable  adjustment  of  .  The 

complete  solution  of  (7-28)  (that  for  (7-29)  is  similar)  is  then  of  the  form: 


[(C)nL,  ..  =  C,  cos  W  t  +  C_  sin  W  t  +  C„  cos  W  t  +  C.  sin  W  t 
(t,i)  1  nj  2  nj  3  n2  4  ng 

(7-34) 

where  C^,  C 2,  C 3  and  are  constants  to  be  determined  by  the  initial  con¬ 
ditions,  and  W  and  W  are  given  by 
’  n«  n„  B  J 


1 

J  e  | 

(- 

J.  4Vt)(U)z[wIE  + (L°",lRs\ 

y  2RS  1 

V  '  s  ) 

(7-35) 


W 


2  R 


(1+  yTT 


4S(int)(i,I)z^WIE 


(Lon)]R 


g 


Let  the  initial  conditions  be  represented  by  ([(C)V]^  .^)o  ,  ([(C)V]^  ^)Q  , 
([(C)N](t  jj)q  ,  ([(C)Nj^.  ^)q  ,  the  y-component  of  the  correction  to  the  indi¬ 
cated  vertical  and  its  first  derivative,  and  the  correction  lo  indicated  north 
and  its  first  derivative,  respectively,  at  the  start  of  the  problem.  The  constants 
CL,  C0,  C,  and  C.  can  be  given  in  terms  of  the  initial  conditions  and  W 

1^0  4  Ilj 

and  W  .  These  constants,  with  (7-34)  and  (7-35),  then  give  the  performance 
n2 

equation.  The  system  is  doubly  periodic,  with  angular  frequencies  W  and 

nl 

W  .  A  similar  result,  with  the  same  initial  conditions,  holds  for  the  y-systeim 
"2 


vlind)y= 


^UntJ(i.i)z 

sec(Lat) 

l(mult)z 

P 

MULTIPLIER 

sin  (Lot) 

MULTIPLIER 
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NOTES 

Hear/  lint  indicates  rigid 
mechanical  connections. 
Light  tint  indicates  signal 
transfer. 


Figure  7-1.  Mathematical  Functional  Diagram  of  Combined  Y-Axis  Vertical  Indi 

cator  and  Z-Axis  Stabilizer. 
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DERIVATION  8 

POSITION  INDICATION  BY  MEANS  OF 
CELESTIAL  LONGITUDE  SIGNAL  COMPOSITION  RATE 


Figure  8-2  shows  the  three-axis  stabilization  system  resulting  from 
combination  of  the  systems  of  Fig.  3-2  in  Derivation  3  and  Fig.  7-1  in  Deri¬ 
vation  7.  Base  motion  isolation  is  also  provided  here  by  the  roll,  pitch  and 
azimuth  drives  at  the  right  of  Fig.  8-2.  The  principal  outputs  of  this  system 
are  latitude  and  longitude,  derived  basically  by  the  resolver  in  the  lower  right 
of  Fig.  8-2,  which  is  changed  from  its  usual  function  to  act  here  as  a  vector - 
component  compositor  (Fig.  -8-1).  Equations  (7-15)  and  (7-21)  in  Derivation  7 
give  the  compositor  inputs: 


Sc  *(int)x  +  *W 


1 

S(cmds)(i,W)x 


1  ^WIE  +  ^  cos  ^Lat^ 


-  [(C)V](t>i)y  [WIE  +  (Lon)]  sin  (Lat) 

+  (lit)  [(ON]  (y)-p[(OV](M)x} 

(8-1) 


Vt.  *  QC0SV> 

RCSOUVCR 

^(out)l 

0 

W0Sin\o.  . 

q(out)2  *  A(fOtr* 

| 

A(rot)  w 

A  1 

M(roO 

II  the  resolver  is  used  os  o  vector  compositor 

%n)l  *  0  cos  A(in)  'houtn  ■  9(ln)l  cos  A(iot)  +  *l(ui)25ln  A(ret) 

q(in)2 ' 0  S,n  A(mt  *>(«««  *  %n>2  cos  A(rot)  “  <lltn)l  s*n  A(rot) 

<Wn  ’  0(cosAmcos  A(tot)  +  sin  A(i„,sin  A(,01))  •  0  cos(A„n>-  A(I<,„)  *  0 


t**"  re  rt 


l  «■ 


'{oum  '  "(rot) 

W°(sin  A((nt  c°s  A(f0tJ  -  cos  A(ltl)sin  A(f01))  ■  0  sin(Ajin)-  A(rotl)  *  °r<c>A(-J 


Figure  8-1.  Resolver  Used  as  a  Vector  Compositor. 
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l  =  i 


z  (malt)z  -  S 


(cmds)(i,W)z 


j-  [WIE  +  (Lon)]  sin  (Lat)  -  (Lat)[(C)V](u)x 


-  C(C)V](t  i)y  (WIE  +  (Lon)]  cos  (Lat)  -  p  t(C)N](y)}  (8-2) 

Under  equilibrium  settled  conditions, 

«C>v)(t,i)x  ■  f<C>VJ(t,i)y=  f(C)Nl(t,t)  ■  0  <8-3> 

The  following  correspondence  can  then  be  set  up  between  the  generalized  vector 
compositor  inputs  of  Fig.  8-1  and  the  compositor  input  signals  of  equations  (8-1) 
and  (8-2). 

A(in)=  A(rot)  =  ^Lat*i  =  A(comp)(rot)(out)  <8'4) 
and  if  sensitivities  are  adjusted  so  that 


S(cmds)(i,W)x  -  S(cmds)(i,W)z 


it  follows  that 


Q  = 

Thus,  the  inputs  are 
1 


(cmds)(i,W)z 


[WIE  +  (Lon)] 


(8-5) 


(8-6) 


q(in)l  S 


(cmds)(i,W)x 


and 


(in>2  '  8(cmds)(i,W)z  IE 


[WIE  +  (Lon)]  cos  (Lat)  =  i^ 


[WT_  +  (Lon)  ]  sin  (Lat)  =  -  i 


(8-7) 


(8-8) 


In  Fig.  8-1,  the  sensitivity  of  the  compositor  for  a  dual  current  input  and 
voltage  output  is  assumed  to  be  unity.  In  the  general  case,  if  this  sensitivity 

is  S(comp)(i  i  -e)’  comPositor  output  is 
1  2 3 


(comp)(i  ,i  ;e)  . 

U  =  e,  w  = - i-=—  WTT,  +  (Lon)  (8-9) 

<out)l  (eompHout)  S(cmds)(iiW)z  IE 
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If  (8-3)  is  not  satisfied,  a  sinusoidal  variation,  constituting  an  inaccuracy  in  in¬ 
dication,  is  superimposed  on  e(comp)(out)  aw*  ^(comp)(rot)(out)'  ^  discussion 

of  the  effects  of  this  condition  will  be  deferred  here  and  taken  up  in  the  next 
report. 


The  component  marked  "clock"  in  Fig.  8-2  furnishes  a  signal  e 


<V 


proportional  to  W^.  The  output  of  the  differential,  which  receives 

« 

e(comp)(out)  anc*  e(W  )  aS  inPuts»  is  proportional  to  (Lon).  Ibis  quantity 
IE 


is  then  integrated.  The  integrator  output  is  proportional  to  the  longitude 
change.  A  longitude  reference  signal  e(Lon)(ref) is  then  added  to  this  to 

obtain  the  angular  difference  between  indicated  longitude  and  the  reference. 
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Figure  P-2.  Mathematical  Functional  Diagram  of  a  Three-Axis  Position  Indicator 
Using  Vector  Composition  of  the  Indicated  Vertical  Angular  Velocity  Components. 


ytW  'f ' 
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DERIVATION  9 


INDICATION  OF  POSITION 
BY  THE  DETECTION  OF  CELESTIAL  LONGITUDE  RATE 
FROM  A  STABLE  PLATFORM 


A.  Latitude  System 


Figure  9-1  shows  a  mathematical  functional  diagram  of  this  system.  In 
the  diagram,  a  single-axis  vertical  indicator  is  shown;  actually,  a  completely 
stabilized  vertical  and  azimuth  indicator  (as  in  Fig.  8-2)  would  be  required, 
with  its  controlled  member  supporting  the  latitude  gimbal  of  Fig.  9-1,  This 
gimbal  supports  the  latitude  gyro  unit,  which,  with  its  accompanying  servo, 
tracks  a  direction  perpendicular  to  the  Earth's  polar  axis.  The  indicated 
polar  axis  will  be  defined  as  (PA).,  the  true  polar  axis  as  (PA)t,  and  the  cor¬ 
rection  to  the  indicated  polar  axis  as: 


t(C)PA] 


(t,i)  "  1(PA)i 


(9-1) 


(See  the  lower  part  of  Fig,  9-1.)  The  torque  summation  equation  for  the  lati¬ 
tude  gyro  unit  is: 

H[[W,E  +  (Lin)][(C)PA)(M).p((C)V](lil)l  sin  (LatJj 

*  sl“  <“>l  +  c°s  <Ut>i]  ■  S(er)<A,M)Ag 

(9-2) 

where  S(er)^  ^) is  t*ie  elastic  restraint  sensitivity  of  the  gyro  unit  and  Ag 

is  the  gyro  unit  gimbal  angle*.  The  microsyn  signal  generator  on  the  gyro 
rotor  shaft  produces  an  output  signal  proportional  to  the  rotor  angle.  This 
signal  is  then  integrated  to  give 


i 


=  S 


1 


(Lat)  "  (sg)(A,e)  (int)(e,i)  p  g 


p(Lat)i  =  S(Lds)(i,A)i(Lat) 


(Laty  o 


(9-3) 

(9-4) 


*  Since  the  dynamic  response  of  the  gyro  unit  is  much  faster  than  the  response  of  the  tracking 
loop,  the  inertia  and  damping  terms  are  neglected. 


When  the  vertical  indicator  is  in  a  settled  equilibrium  state,  so  that, 
approximately, 

[<C>vl(t,i)x  -  f<c)v>(t,i)y  -  t(C)Nl(t,i)  *  0  (9-n> 

the  second  term  on  the  left-hand  side  of  (9-10)  is  zero.  The  solution  of 
(9-10)  is  then,  assuming  a  constant  acceleration  in  latitude: 

r(c)(Ut)](t.)  =  {[<C)(Ut)](M)}  oc„sWnLt  +  1/-{[(CMLat)]M}0si„WnLt 


B.  Longitude  System 

The  torque  summation  equation  for  the  longitude  gyro  unit  is: 

H  {twIE  +  (Lon)t]  cos  [(C)PA](t  i)  +  (p[(C)V](M)x)  cos  (Lat). 

+  (Lat)t  [(C)N](tji)  cos  (Lat).  -  (p[(C)N](t  ,})  sin  (Lat).  +  Wtds 

+  Wlds}  =cpAg  (9-14) 

where  =  -W^.,  the  indicated  angular  velocity  of  the  sidereal  time 
drive  system,  and  Wjds  =  -  (Lon).,  the  angular  velocity  of  the  longitude  gim- 
bal  drive  system  with  respect  to  the  base.  (In  this  analysis  will  be 

assumed  equal  to  WjE),  [(C)PA](t>  ^  and  [(C)N](t  ^  are,  as  before,  assumed 
to  be  small  angles,  nearly  equal  to  their  respective  sines. 

The  gyro  unit  signal  generator  output  is: 


e(sg)  S(sg)(A,e)Ag 


(9-15) 
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From  (9-15)  and  the  definitions  above, 


-  (Lon).  -  Wldg  -  ®(lds)(e>W)  esg 


(9-16) 


which  makes  (9-14): 

(Lon)t  -  (Lon).  +  (p[(C)V)^t  .jx)  cos  (Lat).  +  (Lat)t  [(C)N]^t  ^  cos  (Lat). 

- % -  (Lon). 


-  (p[(C)N]^t  .j)  sin  (Lat).  = 


(lds)(e,W)  (sg)(A,e) 


(9-17) 


Define  [(C)(Lon)j^  ^  =  (Lon)t  -  (Lon). 


where 

and 


(D)(Lon)t  =  (Lon)t  -  (Lon)^Q 
(D)(Lon).  =  (Lon).  -  (Lon)^Q 


(Lon)t  =  (D)(Lon)t  +  (Lon)^ref^t 

(Lon)t  =  (D)(Lon)i  +  (Lon)^e{^  -[(C)(Lon)] 


(ref). 


> (9-18) 


where  [(C)(Lon)]^re^  is  the  correction  to  the  indicated  reference  longitude 
setting.  Therefore, 

p[(C)(Lon)]^t  j  =  (D)(Lon)t  -  (D)(Lon)i  =  (Lo‘n)t  -  (Lon).  (9-19) 
From  (9-17)  and  (9-19), 


P  (1  + 


cp_ 


HS(lds)(e,W)S(sg)(A,e) 


)  [(C)(Lon)],.  ..  = 


(Lon), 


(t,i)  “  HS„  .  W.S,  w.  s  - - 't 

(lds)(e,W)  (sg)(A,e) 


-  (p  [(C)v3(t>i)x>  cos  (Lat)j  -  (Lat)t  [(C)N](t,i)  cos  (Lat)i 
+  (p  [(C)N]^t  .j)  sin  (Lat)j 
When  the  controlled  member  is  in  equilibrium, 


(9-20) 


(9-21) 
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This  reduces  (9-20)  to: 

(p  +  (CT)ldsp2)f(C)(Lon)](ti)  =  (CT)lds(I^n)t  (9-22) 

where 


Ids  HS(ids)(e,W)S(sg)(A,e) 


(9-23) 


The  most  general  solution  of  (9-22)  is 


t 

^icry; 


p[(C)(Lon)]t  =  e  ^CT)ldsj  -j[(Lon)t  e  (CT'lds]  +  [(Lon)^Q  -  (Lon)(.)oJ 

1  (9-24)  J 

Two  cases  are  considered:  the  longitude  rate  may  be  constant  or  changing. 
If  it  is  constant, 


(Lon)t  =  0 


(Lon)t  =  [(Lon)tJo 


(9-25) 


[(C)(Lon)](M)  =  [(C)(Lon)](ref)  +(CT)lds[(Lon)(t)o  -(Lon)(i)o][l  -  e 

(9-26) 

In  (9-26),  [(C)(Lon)]^re^  is  the  correction  to  the  initial  reference  setting  in 

longitude.  The  term  in  (CT)^  is  such  that  a  large  loop  sensitivity  ( (CT)ldg 

small)  results  in  a  highly  damped  system  as  far  as  the  initial  velocity  lag 
([(Lon)tJo  -  [(Lon)jJQ)  is  concerned.  When  (CT)jdg  is  sufficiently  small, 

the  correction  to  the  indicated  longitude  reduces  to 


[(CXLon)]^  =  ((C)(Lon)]^ 


(9-27) 


The  presence  of  an  acceleration  in  longitude  negates  (9-24),  and  the  velocity 
lag  is  then  a  function  of  the  time.  The  lag  may  be  positive  or  negative,  de¬ 
pending  on  the  sign  of  the  acceleration. 


il 


Figure  9-1.  Mathematical  Functional  Diagram  of  a  Position  Indicator  Using  the 
Detection  of  Earth-Rate  from  a  Stable  Platform. 
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DERIVATION  10 

INDICATION  OF  POSITION  BY  COMPARING 
INDICATED  AND  REFERENCE  VERTICALS, 
USING  A  PRE-ALIGNED  INERTIAL  GYRO  UNIT 

A.  Latitude  Indication 

Figure  10-2  shows  a  mathematical  functional  diagram  for  this  system. 
The  function  of  the  controlled  member  drive  system  here  is  the  same  as  in 
Derivation  3,  and  (3-35)  is  applicable. 

Refer  to  Fig.  10-1.  Define 

[(C)(Lat)](t  i)  =  (Lat)t  -  (Lat).  (10-1) 

From  the  figure, 

(Lat).  =  (Lat)t  +  [(C)V](t>.)y  -  [(C)PA](t  .)y  (10-2) 

Therefore 

[(C)(Lat)](tj  #  =  [(C)PAJ(li  ,)y  -  [(C)V](li  ,)y  (10-3) 


witUim  0»  iKHNOiOAT 
IWIUMlSIiliON  UWU10K 

***  **  F.e  H. 


■f +(Lo»lt  +  toPA),titfy  =  -f +(Lo0t+ 
[(C)Lot],tiU5(LoUt-(LoUt 
[(C)Lot)(t  l,= |C)PA)„>Uy  -60 V9(l  ll> 


NOTE  •  This  fiqurt  Is  shorn  os  pi  on  of 
hr  convenience  In  rtprtstntaiion.Ho*trt'f 
this  dots  not  limit  its  opplicobi/ity  to 
tht  etsentioi/y  spherlcoHiigonometric 
pro  bit m  uhdtf  consideration* 


'{<c>v],t(i 


Figure  10-1.  Corrections  to  the  Indicated  Polar  Axis  and  Indicated  Vertical  Pertinent  to 

Latitude  Indication. 


Applying  (3-35)  to  [(C)V]^  ^  in  (10-3) 

[(CKtat)],,.,  =  [(OPA](M)y  -  (t(OV](l.)y,ocosJX  t 

i  s 

<tfC>yJ(t,i)y>0 stn  (10-4) 

s 

Here,  [(C)PA]^.  represents  the  meridian  plane  component  of  a  correction 

to  the  indicated  reference  vertical  due  to  a  possible  misalignment  of  the  in¬ 
ertial  gyro  units  with  respect  to  the  Earth's  true  polar  axis. 

B.  Longitude  Indication 

Refer  to  part  B  of  Derivation  9.  The  case  here  is  similar,  except  for  the 
character  of  the  reference  line  from  which  longitude  difference  is  measured; 
in  the  former  case,  this  line  is  maintained  by  the  combined  action  of  the  con¬ 
trolled  member  stabilization  and  the  orientation  of  the  latitude  gimbal,  while 
in  the  present  case,  the  line  is  maintained  by  the  orientation  of  the  inertial 
gyros  mounted  on  the  inertial  space  gimbal,  and  their  associated  drives. 

Note,  however,  that  the  present  system  requires  an  initial  alignment  of  the 
inertial  gyros,  while  that  of  Derivation  9B  is  self -settling  and  requires  no  such 
alignment. 
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Figure  1 0-2.  Mathematical  Functional  Diagram  of  a  Single-Axis  Position  Indicator 
Using  an  Inertial  Reference  Vertical. 
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